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FOREWORD 

by: Yousry Hamdy 



The St. Marys River sustains a variety of water uses. It is a source 
of water supply to sizeable urban and industrial areas in both Ontario 
and Michigan, is a major tourist attraction and provides a vital 
habitat for many fish species. The Ontario Ministry of the Environment 
is therefore focusing its efforts to permit the greatest number of uses 
based on the best interests of the people of Ontario. 

Several environmental studies, such as those by Hamdy et a1 (1978) and 
Hamdy and LaHaye (1984) have focused on the assessment of the effects 
of Ontario industrial (Algoma Steel Corp. Ltd. and St. Marys Paper), •■ 
and municipal (the easterly and westerly Sault Ste. Marie sewage 
treatment plants) discharges. 

Contaminants such as phenols, amnonia, cyanide, iron and zinc - 
attributable mainly to Algoma Steel discharges, are generally 
restricted to the Ontario shoreline of the river. Levels along the 
Michigan side of the river are similar to those observed at the river 
headwaters. Upon the completion of the redevelopment of Great Lakes 
Power Limited in 1982, the river flow along the Ontario shoreline 
increased from 21% to 36J of the total river flow. This increase is in 
part responsible for the reduction of phenol, ammonia and cyanide in 
the river during recent years. 

Investigations of the river sediments revealed that contaminants such 
as oil and grease and heavy metals are still prevalent along the 
Ontario side of the river. However, data from sediment core samples 
collected during the Upper Great Lakes Connecting Channels Study 
indicates that DDT, PCBs, polynuclear aromatic hydrocarbons, zinc, 
chromium and lead inputs have been reduced over the past 20 years. 

The Upper Great Lakes Connecting Channels Study as well as the St. 
Marys River MISA pilot site investigations will provide further insight 
into the fate and pathways of trace organics in the river ecosystem. 
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Contaminated sediments are a possible source of contaminants to benthic 
organisms and can exert toxic influences on them, either completely 
eliminating benthic populations or reducing the diversity to a few 
tolerant species. However, an Ontario Ministry of the Environment 
study (Persaud et a1 , 1987) found that biological uptake of copper, 
zinc and mercury was evident at locations which exhibited the lowest 
bulk sediment contaminant concentrations. 

The present report on the 1985 benthic invertebrate survey of the St. 
Marys River was prepared by Beak Consultants Ltd. and Is based on 
benthic invertebrate and sediment chemistry data collected from 70 
stations in the river from Point aux Pins Bay to Lake George. 

Characteristics of the benthic community structure of the St. Marys 
River can be described in terms of the following impairment zones: 

1. Severe : 

This zone is found in the Algoma Slip area and in embayments 
downstream from the industrial and municipal discharges along the 
Ontario shoreline of the river. This zone is characterized by 
extreme tubificid dominance. 

2. Moderate : 

This zone, mainly confined to the Ontario shoreline, is 
approximately 500 m wide, extending 4 km downstream from the 
industrial and municipal discharges. Tubificid dominance with high 
densities of nematodes and facultative chironomids ire the major 
characteristics of this zone. 

3. Slight : 

Some recovery was apparent with increased distance from industrial 
and municipal discharges; however, complete recovery was not 
apparent until the lower section of Lake George. Nematode and 
polychaete dominance with moderate densities of tubificids and some 
non-tolerant groups are present. 
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4. Unimpaired : 

This zone was found in the upper reaches of the river and along the 
Michigan side of the river. Communities tended towards chironomid 
dominance, with several non-tolerant groups (e.g., ephemeropterams 
and trichopterans) present, together with low tubi field densities 
and high numbers of taxa. 

It was generally observed that areas with visible oily residues were 
characterized by the absence of ephemeropteran Hexagenia . Substrates 
collected from Lake George in 1985 were contaminated by oily 
substances, similar to the pattern indicated by Hiltunen and Schloesser 
(1983). 

Notwithstanding recent reductions in pollutant loadings, especially _. 
those associated with Algoma Stael Corp. Ltd.. no significant changes 
in the status of benthic communities took place along the Ontario 
shoreline of the river since the earlier studies (1968-1973). 

Algoma Steel Corporation Limited was served with an Amending Control 
Order to limit discharges of phenol, ammonia, cyanide, suspended solids 
and ether solubles (oils and greases). The following outlines effluent 
requirements and associated dates when these requirements should be 
met: 

By March 31, 1990 ether solubles loading must not exceed 1023 kg/d 
and suspended solids will not exceed 5108 kg/d, 

" By June 30, 1989 Algoma will reduce phenol loadings to 22.7 kg/d. 

" By June 30, 1988 ammonia and cyanide levels will be at or below the 
level graphically illustrated by the diagonal line in Schedule 1. 
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SCHEDULE 1 
ALLOWABLE CONCENTRATIONS OF TOTAL AMMONIA AND FREE CYANIDE AT THE TERMINAL BASIN OUTFALL 
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SUMMARY AND RECOMMENDATIONS 

A study was undertaken in 1985 to examine spatial and temporal trends in the structure 
of the benthic community of the St. Marys River in relation to sediment contamination 
and wastewater sources in the vicinity of Sault Ste. Marie, Ontario and Sault Ste. Marie, 
Michigan. This survey concentrated on the lower river (Lake George, Little Lake 
George, Lake Nicolet) and in the upper river in the vicinity of Algoma Steel, as mid-river 
locations had been relatively intensely sampled in 1983. Some of the stations sampled in 
1983 were resampled in 1985 to determine the compatibility of the two data sets. 

The water quality of the St. Marys River has been impaired along the Canadian shore by 
discharges of organic and inorganic contaminants from pulp and paper and steel mill 
industries (St. Marys Paper and Algoma Steel, respectively), and the municiparsewage 
treatment plant (STP). Benthic surveys in 1968 and 1973 demonstrated zones of severe 
impairment along the Canadian shoreline downstream of the various discharges. 
Subsequent to these surveys, considerable progress has been made by industry in reducing 
contaminant loadings to the river. The 1983 benthic survey indicated that while some 
minor improvements had occurred, zones of severe impairment were still apparent along 
the Canadian shoreline. 

The benthic survey detailed in this report was carried out in September/October 1985. 
Surficial sediment samples were collected concurrently and submitted for physical and 
chemical analyses, and a summary of the analytical results is provided in this report. A 
detailed listing of the results of these analyses will be reported in a forthcoming 
document prepared by the Ministry of the Environment. Benthic samples were collected 
in triplicate at 70 stations between Point aux Pins Bay and Lake George. Most stations 
were sampled using a Ponar grab and a 200 um mesh sieve, although, at a few stations, 
the sediments were too coarse and were sampled by an airlift sampler/200 um mesh sieve 
combination. Additional samples were collected at the 11 stations common to the 1983 
and 1985 surveys to evaluate the effects of sampling method (ponar or airlift) and mesh 
size (200 um vs 500 um mesh) on the results. Benthic organisms from each sample were 
sorted from the debris retained by the sieve, identified and enumerated. 
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Data from the 1985 St. Marys benthic survey were not amenable to pooling with those of 
the 1983 survey, although the surveys were conducted at the same time of year. Total 
densities were found to be generally greater in 1985, more so at some stations than 
others. Different sorters and sorting methods were identified as significant contributing 
factors. A principal components analysis suggested that a faunal shift in species 
associations, either real or sorter-related, had also occurred between the two survey 
years. 

As in 1983, the tubificid community associated with outfall areas was still apparent in 
1985. Nematodes also reached high densities in all but the most severely enriched 
locations. General reductions in pollutant loadings from the Algoma Steel terminal basin 
and the St. Marys Paper operations since 1973 appear to have contributed to-minor 
community changes, but were not reflected in major improvements in the pollutional 
status of the benthic community. 

Cluster analysis based on species composition in 1985 showed that benthic communities 
characterized by tolerant species occurred in areas downstream of discharges, and other 
communities characterized by intolerant species occurred in areas remote from point 
source effects. Similar impact zones can be derived from a community index reflecting 
dominance of nematodes and immature tubificids without setae. Discriminant analysis 
identified a heavy metal-particle size gradient and a pesticide-particle size gradient m 
sediment quality, which together provide a basis for separation of impacted station 
clusters from each other and from unimpaired stations. 

Impacted stations generally occurred downstream of the Aigoma Steei, St. Marys Paper 
and easterly Sault Ste. Marie sewage treatment plant (STP) discharges. Some recovery 
was apparent with increasing distance downstream from the Algoma Steei and St. Marys 
Paper discharges; however, complete recovery was not apparent until the lower section 
of Lake George. The combined impact of effluents from all sources was still evident in 
Lake George, some 2k km downstream from the Algoma Steel and St. Marys Paper 
discharges. In general, impacted locations were restricted to the Canadian portion of the 
St. Marys River, with a clean water fauna characterizing the relatively unindustrialized 
Li.S. shore and all portions of the river upstream of pollution sources. 
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The distribution of Hexagenia nymphs was found to be influenced by the occurrence of 
visible oil in the substrate of the St. Marys River, Lake Nicolet, and the Lake George 
Channel downstream into Lake George. Generally, absence or low densities of nymphs 
coincided with the presence of oil in the sediments. The pattern observed in 1985 for the 
distribution of oil and HexaRcnia was similar to that reported in 1975, which indicated 
the occurrence of oil as far downstream as the upper half of Lake George. 

Based on the results of this survey, three recommendations for future environmental 
monitoring can be made: 

1. A benthic survey on the St. Marys River should be repeated in 
approximately 5-10 years to document any changes in the benthic 
community in response to on-going pollution abatement programs. 

2. In view of the limited benthic community response to pollution abatement 
over the past few years, analysis of temporal trends in sediment quality, 
and of direct organism responses (e.g., bioaccumulation, toxicity), should be 
undertaken. (Such studies were initiated in 1986.) 

3. In order to improve data comparability between benthic surveys a set of 
basic guidelines regarding sampling design, sampling methods and sorting 
techniques should be developed. In addition, QA/QC requirements 
applicable to benthic surveys should be developed to check for systematic 
bias between laboratories and sorters. This recommendation applies to 
benthic surveys in general. 
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SOMMAIRE ET RECOMMANDATIONS 

En 1985, on a entrepris une etude dans le but d'ezaniner, dans 
1 ' espace et dans le temps, les tendances relatives a la 
composition du benthos dans la riviere St . Marys au regard de la 
contamination par les sediments et des sources d'eauz usees dans 
les environs de Sault Ste. Marie (Ontario) et de Sault Ste. Marie 
(Michigan). Cette etude a porte principalement sur la basse 
St. Marys (lac George. Petit lac George, lac Nicolet) et sur la 
haute St. Marys dans les environs de la societe Algoma Steel, 
etant donne que le cours moyen da la riviere avait fait I'objet 
de prelevements relativement nombreux en 1983. On a effectue en 
1985 des prelevements a des postes qui avaient deja ete etudies 
en 1983 afin de determiner la compatibilite des deuz series de 
donnees obtenues. 

La qualite de 1 ' eau de la rlviSre St. Marys s'est deterioree le 
long des rives canadiennes en raison de deversements de polluants 
organiques et inorganiques provenant de I'usine de pates et 
papiers St. Marys Paper, de 1 ' acierle Algoma Steel et de I'usine 
d'epuration des eauz d'egout de Sault Ste. Marie. Des etudes 
benthiques menees en 1968 et 1973 ont revele la presence de zones 
fortement polluees le long des rives canadiennes en aval du lieu 
des deversements. A la suite de ces etudes, les pollueurs en 
cause ont pu reduire considerablement leurs deversements dans la 
riviere. L' etude benthique de 1983 a revele que malgre certaines 
ameliorations mineures, des zones fortement polluees subsistaient 
le long des rives canadiennes. 

L' etude benthicjue dont les details figurent dans le present 
rapport a ete effectuee en septembre et octobre 1985. Des 
echantillons de sediments superficiels ont ete preleves et soumis 
k des analyses physiques et chimiques dont les resultats sont 
resumes dans le present rapport. La lists detaillee de ces 
resultats figurera dans un document ulterieur du ministere de 



1 ' Environnement . Des echantillons de benthos ont ete preleves h. 
trois reprises a 70 postes entre la bale Pointe auz Pins et le 
lac George. A la plupart des postes, les prelevements ont ete 
effectues a I'aide d'une pince Ponar et d'un tamis a mailles de 
200 um: a certains postes, les sediments etaient trop grossiers 
et il a fallu utiliser une suceuse h. air combiuee a un tamis a 
mailles de 200 um. On a preleve des echantillons supplementaires 
aux II postes communs a 1' etude de 1983 et a celle de 1985 afin 
d'evaluer les effets de la metliode d' echantillonnage (par pince 
Ponar ou suceuse a air) et du tamis (200 urn ou 500 um) sur les 
resultats. Les organismes benthiques de ctiaque echantillon 
etaient extraits des debris retenus par le tamis. puis identifies. 
et denombxes. 

Les donnees de 1' etude bentiique de la St. Marys effectuee en 
1985 n'ont pu etre combinees a celles de 1' etude de 1983, bien 
que les deux etudes aient ete menees environ a la meme periode de 
I'annee. Les densites totales. en particulier h certains postes, 
se sont revelees plus elevees en 1985. A cet egard, 1 ' utilisation 
de materiel d' echantillonnage et de methodes de triage differents 
ont ete consideres comme des facteurs determinants. L' analyse 
des composantes principales a permis de constater qu'une^ 
modification des associations d'especes, reelle ou causee par le 
triage, s'etait produite entre les deuz etudes. 

Comme en 1983, la presence de Tubiflcidae dans les zones de 
deversement etait decelable en 1985. Les nematodes se trouvaient 
egalement en forte densite. sauf dans les zones les plus 
polluees. La reduction generale de la quantite de polluants 
deverses du bassin terminal de I'acierie Algoma Steel et de 
I'usine de la societe St. Marys Paper depuis 1973 semble avoir 
entraine certains changements mineurs, mais non pas de diminution 
Importante de la pollution du benthos. 

En 1985. une analyse typologique fondee sur la composition des 
especes a montre que les communautes benthiques comportant des 



especes resistaixtea a la pollution se trouvaleat dans lea zones 
aituees en aval des points de deversement , et que les cominunautes 
composees d' especes non reslstantes se situaient dans les zones 
eloignees des sources ponctuelles. Des zones toucliees semblables 
peuvent etre determinees a partir d'un indice refletant la 
preponderance des nematodes et des Tubificldae sous-developpes 
sans soles. L' analyse discriminante a permis de caracteriser les 
sediments selon la taille des particules de metaux lourds et de 
pesticides cju'ils contiennent, ce qui constitue un moyen de 
distinguer les postes touclies h la foia ©ntre eux et par rapport 
aux postes non touches. 

La plupart des postes touclies etaient situes en aval des points 
de " deversement de I'acierie Algoma Steel. de I'usine 
St. Marys Paper et de I'usine de traitement des eaux d'egout de 
Sault Ste. Marie, h. I'est. Plus loin en aval de I'acierie et de 
I'usine de pates et papiers. une certaine diminution de la 
pollution a ete constatee; cependant . ce n'est qu'a partir de la 
section sud du lac George que la qualite de 1 ' eau revenait k la 
normale. L' incidence combinee des effluents de toutes sources 
§tait encore decelable dans le lac George, a environ 24 km en 
aval des points de deversement de I'acierie et de I'usine de 
pates et papiers. Dana 1' ensemble, les zones touchees etaienfr 
situees dans la partie canadienne de la riviere St. Marys (une 
faune d'eau non polluee se trouvait le long des rives 
etatsuniennes. relativement peu industrialisees) et dans toutes 
les parties de la riviere en amont des sources de pollution. 

La repartition des nympiies d ' Hexagenla s'est revelee influencee 
par la presence d'huile visible dans le substrat de la riviere 
St. Marys, du lac Nicolet et du canal du lac George en aval 
jusqu'au lac George. L' absence ou la faible densite de nymphes 
correspondait generalement a la presence d'huile dans les 
sediments. La repartition de I'huile et des Hexagenia observee 
en 1985 etait semblable a celle qui avait ete relevee lors d'une 
etude effectuee en 1975 et qui avait etabli la presence d'huile 
en aval Jusque dans la premiere moitie du lac George. 



Des resultata de 1' etude. decoulent trois recommandations 
relatives a la surveillance environnementale : 

,1. II convient de refaire une etude benthique de la 
riviere St. Marys dans 5 a 10 ana afin de determiner si 
les programmes d' assainissement en cours ont entraine 
une modification du bentlios. 

2. Etant donne le peu d' amelioration du benthos constate 
depuis quelques annees malgre les mesures 
d 'assainissement, il convient d'entreprendre 1' analyse 
des tendances temporelles de la qualite des sediments 
et 1' analyse des effets directs sur les organismes 
(p. ex., bioaccumulation, tozicite). Des etudes de ce 
genre ont ete entreprises en 1986. 

3. Afin d'auneliorer la comparabilite des donnees des 
differentes etudes benthiques. il convient d'elaborer 
un ensemble de lignes directrices f ondamentales portant 
sur le plan et les methodes d' echantillonnage et sur 
les tecliniques de triage. En outre. il y a lieu 
d'etablir des normes d' assurance et de controle de la 
qualite touchant les etudes benthiques afin de verifier 
s ' il se produit des hearts systematiques entre les 
laboratoires et les appareils d' echantillonnage. Cette 
recommandation s' applique aux etudes benthiques en 
general. 



1.0 INTRODUCTION 

1.1 Background and Objectives 

The maintenance of a high level of water quality in the St. Marys River is one of the 
major concerns in the upper Great Lakes system. The river is a source of water supply 
for a sizeable urban and industrial area in both Ontario and Michigan, is a major tourist 
attraction and provides a vital habitat for many important fish species. Therefore, 
degradation of water quality can have serious implications for water use in both Canada 
and the United States. 

Pollution resulting from steel and paper mill discharges has been documented since the 
19'fO's and has been manifested in impairment of water quality, sediments and the 
zoobenthos. The Ontario Water Resources Commission (OWRC), in cooperation with the 
International 3oint Commission (DC), conducted a survey of the impacted benthic 
community in 1967 (Veal, 196S). Considerable progress has, however, been made since 
1970 by Algoma Steel Corporation in reducing phenol discharges, by Abitibi Paper (nov. 
St. Marys Paper) Company Ltd. in reducing suspended solids loadings, and by Sault 5te. 
Marie, Ontario in improving the removal of settleable organic matter from the municipal 
sewage treatment plant. 

As part of Ontario's contribution to the three-year International Reference Study on the 
pollution problems of Lake Huron and Lake Superior, the Ontario Ministry of 
Environment (MOE) intensified assessment activity in the river in 1973 and 197^* for the 
purposes of: 

determining the water quality status in light of improvements In waste 

treatment; 
predicting the ultim.ate effectiveness of ongoing remedial measures in 

restoring river water quality; and 

investigating the need for further remedial actions. 

Information was also required to assess water quality in terms of DC Agreement 
objectives and to investigate the transboundary movement of pollutants. 



Results of a 1973 benthic survey (Hamdy _et _aL, 1978) indicated little improvement since 
1967 in the benthic community. Zones of severe or strong impairment were identified 
adjacent to and downstream of the Algoma Steel, Abitibi-Price and Sault Ste. Marie STP 
discharges (Figures Al.l and A 1.2, Appendix 1). As a result of this study, further 
reductions in waste loading and additional follow-up studies were recommended. 

lEC Beak Consultants Ltd. (lEC BEAK), under contract to the MOE Great Lakes Section 
undertook a benthic and sediment survey of the St. Marys River in the autumn of 1983 
(McKee et. aL, 198'^). The purpose of this survey was to examine sediment 
quality /ben thological relationships in light of changing hydrological conditions brought 
about by recent reconstruction of the Great Lakes Power Corporation generating station 
and ongoing pollution abatement programs at Algoma Steel. 

It was concluded that the benthic communities of the St. Marys River m 19S3 were 
generally similar to those found in 1968 and 1973, with zones of severe or moderate 
impact adjacent to and downstream of the pulp and paper mill, steel mill and sewage 
treatment plant discharges. A severely impacted community characterized by high 
densities of Limnodrilus/Tubifex reported earlier was less pronounced in 1983, and the 
severe impact zone was confined to a region downstream of St. Marys Falls (Figure A1.3, 
Appendix 1). This has been identified as a region of heavy metal contamination. 
Reductions in pollutant loadings by industry appeared to have contributed to the minor 
community changes. No impacts from transboundary transport of contaminants were 
apparent in the benthic community along the Michigan shoreline. 

The 1983 survey was less extensive than in earlier years, with little or no sampling 
downstream of the industrial region in Lake George, Little Lake George or Lake Nicolet, 
or upstream of Algoma Steel to Vlhitefish Bay. Since benthic community effects were 
still apparent throughout the industrial region, additional sampling in the upstream and 
downstream areas was recommended to define their full extent. On the basis of a survey 
by Hiltunen and Schloesser (1983) in 197'*-75, the benthic community is believed to be 
impacted to a point at least 30 km downstream of Sault Ste. Mane due to oiiy substances 
in the sediments. 

Accordingly, the Ministry of Environment (Water Resources Branch, Great Lakes Section) 
sponsored a 1985 survey to augment the 1983 results. While emphasis was placed on the 
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lower and upper rivers, some of the 1983 stations were resampled in 1985 to determine 
the comparability of 1983 and 1985 data. This report documents the results of the 1985 
benthic survey, compares the 1985 data with those from the earlier studies, and 
evaluates relationships between sediment quality and benthic community structure. 

1.2 Description of Study Area 

The St. Marys River is the connecting waterway between Lake Superior and Lake 
Huron. The upper portion of the river extends from Whitefish Bay to the St. Marys 
Falls. A series of works (locks, power diversion channels, regulating gates) have been 
constructed as aids to control the outflow of Lake Superior for navigation and for 'power 
development. 

The lower river has an irregular shoreline and contains three large islands, Sugar and 
Neebish Islands on the U.S. side and St. Joseph Island on the Canadian side. Sugar Island 
diverts the river below the locks into two channels, Lake Nicolet and Lake George 
(Figure 2.1A). St. Joseph Island subdivides the flow from Lake George into the middle 
and east channels, while Neebish Island subdivides the flow from Lake Nicoiet into 
middle and west channels. 

Geology and Topography 

The rock formations of the land bordering the St. Marys River consist of schists, granite 
and diabases of Precambrian age which are mantled and in places deeply buried 
underneath glacial drift. The area on the Canadian side which borders the St. Marys 
River is a drift-covered strip from one to five miles in width. The land surface rises 
toward the north in a succession of terraces and reaches a maximum altitude of about 12C 
m above the mean water surface elevations of both Lake Superior and Lake Huron which 
are ISO m and 175 m above sea level, respectively. On the United States side, the 
topography is characterized by gradually rising terrain which develops into a flat area 
with an altitude of about 90 m above the bordering lakes. 

Water Uses 

Domestic later Use - The Upper St. Marys River is the source of supply for the U.S. and 
Canadian cities of Sault Ste. Marie. The total average pumpage from the river is about 
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O.t m^/s, from a new intake situated at Gros Cap, Ontario and about 0.2 m /s for Sault 
Ste. Marie, Michigan from the river upstream of the locks. No public water supply is 
taken from the lower St. Marys River. The majority (approximately 70%) of Sault Ste. 
Marie, Ontario is supplied by wells. 

Indtestrial Water Use - Algoma Steel and the St. Marys Paper Company withdraw about 
6.'* and 0.8 m^/s respectively from the river. 

Hydroeiectric Power - The Great Lakes Power Corporation Plant (Ontario), located 
between the Canadian lock and the mainland, used approximately 400 to 500 m /s in 1973 
and 930 m-'/s in 1983, while the hydroelectric plant located between the U.S. locks and 
the rapids requires about 350 m^/s. The Edison Sault (Michigan) hydroelectric; power 
plant is served by a «f km long canal which diverts water from a point just above the U.S. 
locks and delivers it to the plant. 

Navigation - Vessel traffic through the locks has been heavy. However, the average 
number of vessels passing through the locks has decreased from 26,122 vessels in 1953 to 
12,712 in 1970, 8,510 in 1985 and 8,3<t5 in 1986. In recent years, shipping tonnage has 
fallen from the 1970 level of about 86 million tons of cargo to 7«* million tons in 1985 and 
70 million tons in 1986 (Corps of Engineers, pers. comm., June 1987). The vessels carry 
mainly crude oil, grain, steel, coal and petroleum products, taconite and iron ore between 
Lake Superior and the industrial centres on the lower lakes. 

Recreation - The bays and lakes in the lower river provide areas with excellent 
recreational possibilities. Park use is intensive in the two parks nearest the cities of 
Sault Ste. Marie. There is private bathing along all parts of the river other than the 
immediate vicinity of the cities. The Algoma Health Unit restricts swimming in the 
river from the locks to Bells Point at Little Lake George. Boating and fishing are 
popular forms of recreation throughout the entire area. 

Waste Disposal 

Locations of municipal and industrial discharges into the river are shown in Figure 2.1 A. 
Municipal waste waters from both cities are treated in primary plants with chlorination 
facilities. The effluent from the City of Sault Ste. Marie, Ontario is discharged to the 
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Lake George Channel at a rate of about C.6 m-^/s, from the existing eastern faciht>. A 
new westerly plant (which came into operation in early 1986) discharges to the river 
upstream from the Algoma slag dump site at a rate of O-Z** m'^/s. The U.S. facility is the 
ajor U.S. point source discharger to the international section of the river, at a rate of 



m 



0.13 m^/s. 

Industrial waste discharges from Sault Ste. Marie, Ontario total about 7.2 m /s, about 
90% of this from Algoma Steel. There are no industrial discharges from Sault Ste. Mane, 
Michigan to the river. 

Trend analyses for the effluent qualit> of the Algoma Steel terminal basin (which 
discharges about 80% of Algoma Steel flows) indicated that effluent concentrations from 
Algoma* were reduced substantially between 1972 and 1985 (Figures Al.i^ tcJ'Al.S, 
Appendix 1), as follows: phenolics (80%', ammonia (80%), cyanide (98%) and sulphides 
(99%). Total suspended solids loadings increased slightly (10%) in 1985 over 1972-73 
levels, after decreasing substantially in 1982-83 from levels recorded in 1974-81. These 
reductions resulted in concentrations in the St. Marys River approaching the 1978 Great 
Lakes Water Quality Agreement (GLVIQA) and Provincial Objectives (PUQ) for the 
protection of aquatic life (MCE, unpublished data). This decrease was likely attributable 
to the reduction in point source loadings, as well as an increase (approximately two-fold) 
in river flow along the Ontario shoreline resulting from increased diversion to Great 
Lakes Power in 1982. In 1985, flow rates from Lake Superior were the highest on record 
at 3,766 m^/s (as compared to a mean flow of 2,1*0 m^/s recorded over 1955 to 1979), 
with about *0% of this diverted from Lake Superior through the Great Lakes Power 
facility. These higher flows would tend to increase dilution of the discharges as well as 
to displace existing contaminated bedloac sediments. 

The other main industrial waste water discharge is from the St. Marys paper mill 
operation. This discharge takes place through the main mill sewer and the sulphite 
sewer. The latter, however, has been terminated. The total flow from the main mill 



G. LaHaye, MOE, pers. comm. 
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sewer was 0.8 m^/s tn 1973. Effluent loadings of suspended solids and BODj from St. 
Marys paper mill have been declining steadily since I96S. Suspended solid loadings were 
reduced 75% between 1973 and 1983. 



2.0 METHODS 

2.1 Field Survey Methods 

The 1985 survey was carried out from 2U September to 04 October 1985 using the 
Ministry of the Environment survey vessel, "Monitor V". Stations were carefully located 
using the vessel's radar positioning equipment. Sampling locations are illustrated in 
Figure 2.1A and 2. IB. A total of 71 stations were sampled in 1985, with 11 stations 
corresponding to 1983 locations. 

The station numbering system used started from stations sampled at the south end of 
Lake George, increasing upstream. Where transects were sampled, the number system 
denoted the distance upstream or downstream from the St. Marys Falls and the distance 
from the U.S. shore. For example, Station SMD 6.fE-'t90 indicates a location 6.!> mi 
downstream of the falls in the eastern channel and 490 m from the U.S. shore. 

A preliminary survey completed in 1983 (McKee et aL, 1984) indicated that three 
replicate samples at each sampling station were sufficient for estimation of log- 
transformed total organism density with a standard error less than or equal to 20% of the 
mean density. This requirement was not likely to have changed since 1983. Therefore, 
three replicate samples were collected at each 1985 sampling station. 

Benthic samples were collected primarily by Ponar, with two stations requiring airlift 

sampling due to substrate characteristics. The Ponar grab sampled a bottom area of 0.05 

m^. The diver-operated airlift sampler was equipped with a 6 cm diameter opening at 

2 
the intake, and 200 urn mesh Nitex collecting bags. A 0.075 m circular quadrat was 

used on the river bed to demarcate the area sampled. Invertebrates were lightly brushed 

by hand by diver from any rock and cobble in the sampled area and drawn into the intake 

of the airlift. All benthic samples were washed through 200 urn mesh and were preserved 

in 5 to 10% formalin. 

Benthic samples were then dispatched to the laboratory for sorting, identification and 
enumeration of benthic fauna. Benthic organisms were sorted from sediment and debris 
using a dissecting microscope (lOx power) and grouped by major taxa. After sorting, a 
subsample of the tubificids (> 20%) and chironomids (> 10%) in each sample was cleared 
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and mounted in "CMCP-10" mounting medium on microscope slides. Chironomids were 
generally decapitated to facilitate clearing of the head capsule. Benthic organisms were 
identified using various taxonomic keys - Klemm (1985) for Oligochaeta; Oliver and 
Roussel (1983), for Chironomidae; Mackie et aL (1980) and Clarke (1981) for MoHusca; 
Wiggins (1977) for Trichoptera; and Pennak (1978) for other groups. All abundance data 
were expressed on a 0.05 m' basis. 

Surficial sediments were collected at each station for physical/chemical analysis using a 
Shipek sampler. This device samples a 0.0'f m area of bottom and was subsampled to a 
depth of 3 cm. Depending on substrate type, three or more sediment samples were 
needed to provide sufficient sample volume to meet the analytical requirements. MOE 
also meaured Eh (redox potential) of the samples before they were removed from the 
Shipek sampler. The presence of any odour (i.e., phenolic) or oily sheen was noted, and 
sediment texture and colour were described. The presence of aquatic plants and woody 
debris was also recorded, along with water depth and perceived current. Sediments from 
each station were homogenized in a clean, solvent-rinsed stainless steel tray using 
stainless steel utensils. Sediments were then partitioned into subsamples for analysis of 
trace organic and inorganic constituents by the MOE laboratories in Rexdale. Results of 
the sediment analyses will be presented in a forthcoming report from the Great Lakes 
Section, Water Resources Branch. A summary of the data is presented herein for the 
purposes of examination of sediment quality-benthic community structure relationships. 

2.2 Statistical Analysis 

2.2.1 Comparability of 1983 and 1985 Benthic Surveys 

A number of preliminary analyses were performed to assess the comparability of I9S3 
and 1985 survey data. Comparability would permit pooling of 1983 and 1985 data to 
provide a larger database for further analysis. Potential sources of data set 
incompatibility include: 

1. differences in sampling equipment. - most stations in 1983 (35 of hZ 
surveyed) were sampled by airlift sampler, while the others were sampled 
using a Ponar grab. In 19S5, all but two stations were sampled by Ponar. A 
200 um mesh was used for sieving in both years; 
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2. differences in sorting efficiency - changes in personnel between surveys 
could produce differences in retrieval of organisms from 1983 and 1985 

samples; and 

3. real changes in benthic faunal densities and associations between 1983 and 

1985. 

Each of these questions was separately addressed. The results are summarized below, 
with supporting data and statistics in Appendix 2. 

Sampling Equipment 

Three stations were sampled by multiple methods during the 1985 survey (SMD 5.0E-6C, 
SMU 1.0-2130 and SMU 2.5-1830). The methods included airlift sampling with a 200 um, 
mesh, Ponar sampling with a 200 urn mesh, and Ponar sampling with a 500 um mesh (a 
method used in St. Marys River surveys prior to 1983). Five replicate samples were 
collected at each station by each method, and total organism densities were compared 
(Table A2.IA). 

Two-way factorial ANOVA on log-transformed densities indicated a significant 
interaction between station and method effects (p = 0.05) when Ponar 500 um mesh data 
were included (Table A2.1B). This interaction disappeared, and methods did not differ 
significantly, when Ponar 500 um mesh data were excluded (Table A2.1C). Thus, Ponar 
and airlift sampling methods should agree as long as the same mesh size (200 um) is used 
in sample sieving. 

Sorting Efficiency 

Four samples collected in 19S5 (118, SMD5.0E-210, -370 and -520), three of which were 
also sampled in 1983 (SMD5.0E-210, -370 and -520), were independently sorted (at the' 
subsample level) by two different individuals in 1985. One sorter attempted to reproduce 
the sorting method used in 1983 (i.e., without the aid of a dissecting microscope), whiie 
the other sorter used a dissecting microscope at about lOx power (as used for all other 
samples sorted in 1985). Total organism densities of the subsample portions reported by 
these sorters, and by the 1983 sorter, for each of three replicates per station were 
compared (Table A2.2A), and are presented in Figure A2.1. Two-way factorial ANOVA 
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on 1985 log- transformed densities indicated a significant interaction between station and 
sorter (Table A2.2B). Thus, sorter effects vary from one station to another. They could 
represent a source of database incompatibility, and it would not be possible to adjust one 
year's data to be comparable with those of another year. Sorting done with the aid of a 
microscope usually resulted in a greater recovery of organisms than sorting done without. 

Changes in Benthic Density 

The possibility of changes in total organism density between surveys was assessed at 11 
stations sampled both in 1983 and 1985 (SMD 5.0-60, -210, -370, -520, -640; SMU 0.5-200, 
-300; SMU 2.5-1300, -1830, -2£j70; and SMU 1.0-2130). Detailed results from these 
locations are presented in Table A2.3A (Appendix 2). Two-way factorial ANOVA on log- 
transformed densities indicated a significant interaction between station and year 
(p < 0.05; Table A2.3B). Densities were always higher in 1985, but the magnitude of the 
difference varied strongly from one station to another. 

Density increases between surveys could represent real faunal changes or they could be 
related to the use of different sorters in different surveys. In either case, because the 
effect varies among stations, it is not possible to adjust one year's data to be comparable 
with those of another year or to combine the two data sets. 

Changes in Faunal Association 

Species associations were compared between surveys at the 11 stations sampled in 1983 
and 1985. A Principal Components Analysis (PCA) was performed on each year's data 
using both log-transformed and untransformed data. This analysis is based on the matrix 
of correlations between densities of different species. The most distinctive group of 
associated species is identified as the first principal component of the benthic 
community. Secondary species groups, distributed independently of the others, are 
identified as subsequent principal components. 

Patterns of species association, as represented by the plot of species 'loadings' on the 
first two principal components, were quite different between surveys, whether 
transformed or untransformed data were used (Table A 2.4 and Figure A 2.2 present 1983 
results and Table A2.5 and Figure A2.2 present 1985 results, based on untransformed 
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data). Moreover, only about one-third of the key species contributors to the first 
principal component were common to the 1983 and 1985 surveys. This difference 
between surveys suggests that real or sorter-related faunal shifts have occurred, and that 
1985 data should not be pooled with 1983 data for analysis of the benthic community. 

Immature tubificJds were treated as two distinct taxa (capilliform and non-capilliform) in 
principal components analysis, since at most stations where they occurred in 1985, 
mature tubificids were absent, preventing any logical representation of tubificid species 
composition. This convention was applied to both the 1983 and 1985 data sets for PCA. 

2.2.2 Cluster Analysis of Spatial Patterns 

The 1985 benthic communities were defined by means of cluster analysis (Anderberg, 
1973; Green, 1979; Gordon, 1981). This family of techniques for pattern recognition 
identifies stations with similar benthic species composition, and groups biologically 
similar stations into clusters. The species characterizing each cluster of stations are 
interpreted as distinctive natural communities. In cases where pollution is an important 
determinant of community structure, these communities tend to be spatially distributed 
in relation to pollution sources. 

Several cluster analysis techniques were used in order to identify robust spatial patterns 
which are evident by several techniques. The techniques used were as follows: 

1. convergent K-means partitioning, 

2. sum of squares agglomeration (Ward's method), 

3. group average link agglomeration, and 
k. Howard-Harris polythetic division. 

The first technique is non-hierarchical and assumes a priori knowledge of the number of 
station clusters (i.e., faunal assemblages) in the community. As this is seldom known, a 
range of values for the number of clusters was tried, and the solutions were compared. 
Methods 1, 2 and 3 were performed using the SPSS statistical package for the IBM-XT 
(Norusis, 1986). 

Methods 2 and 3 are hierarchical agglomerative techniques in which stations or station 
groups are sequentially fused according to their similarities. The process continues until 



- 16 



all stations are accounted for, producing a solution (i.e., station grouping pattern) at each 
fusion stage. Both methods are based on inverse measures of station similarity, either 
Euclidean distance (D) or Squared Euclidean distance (D ). The distance formula for 
Stations j and k (Norusis, 1986) is: 



I I 

d2 (j, k) = ^ (X^j - Xhk)2 

where: X^j = log-transformed abundance of species h at Station j. Method 2 (Sard's 
method) was based on D and Method 3 on D. 

Method ti (Harris, 1985) is a hierarchical divisive technique, using the first fVincipa! 
component of species composition as a splitting criterion to sequentially divide the 
stations into increasing numbers of groups. At each stage, stations are grouped so as to 
minimize within-cluster variance of the principal component (i.e., maximize within- 
cluster similarity). Statistical software was provided by B. Harris, University of 
Pennsylvania. 

The cluster solution producing the greatest degree of community structure was chosen to 
represent the St. Marys River benthic community. Structure is indicated by well-defined 
station clusters, internally homogeneous and dissimilar from each other. Solutions of this 
nature tend to be the most interpretable. Figure A 3.1 (Appendix 3) illustrates the 
relationship between number of clusters and maximum biological distance between 
clusters for the hierarchical agglomerative techniques, and demonstrates the superiority 
of the Ward's solution. This solution shows the greatest instantaneous rate of increase in 
between-cluster distance, and suggests a break point or sharp increase in communit> 
structure at the seven-cluster level. Thus, the seven-cluster Ward's solution was chosen 
to represent spatial patterns in the benthos. 

Other cluster analysts solutions are compared to the Ward's solution in Table A3. 1 
(Appendix 3). This comparison shows a good concordance between the seven-cluster 
Ward's solution and the seven-cluster Howard-Harris solution (52 of 70 stations falling 
within the same station clusters and, following cluster characterization, agreement on 
impacted vs. unimpacted status of 66 stations). Thus, the representation of 1985 spatial 
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pattern appears to be robust, and the cluster analysis methods used in 1983 and 1985 
agree closely (see McKee et aL, [3U}. 

Immature tubificids were treated as two distinct taxa (capilliform and non-capilliform) in 
cluster analysis of the 1985 data, since at most stations where they occurred, mature 
tubificids were absent, preventing any accurate representation of tubificid species 
composition. 

This approach differed from the approach used in 1983 when immature tubificids were 
usually accompanied by mature individuals, and only the mature individuals were included 
in the data set. 

2.2.3 Temporal Comparison of Spatial Patterns 

Characterization of station clusters in 1985 was based on comparison of total organism 
density, species richness, tubificid dominance and species composition between clusters. 
This information was used to identify impacted station clusters, and from the spatial 
distribution of these clusters to delimit impact zones. Spatial patterns and impact zones 
in 19S5 were related to known past and present pollution inputs, and compared to spatial 
patterns and impact zones during previous surveys (McKee et^aL, 198^), 

I.IA Sediment Quality Relationships of Benthic Communities 

Several independent methods were used to examine the sediment quality characteristics 
of the St. Marys River benthic communities. First, discriminant analysis was performed 
on the sediment quality data in order to identify environmental gradients in sediment 
quality which tend to distinguish between the stations in different benthic fauna! clusters 
(Section 2.2.2). This method has been described by Green (1979), and other authors (e.g., 
Hutchinson, 1978), as a means of representing the ecological niche space occupied by 
different communities. 

A second method of examining the sediment quality - benthic faunal relationships is by 
derivation of a quantitative benthic community index, with a value at each station, to be 
used as a dependent variable in regression analysis of the sediment data. Thus, the 
sediment characteristics which are the best predictors of the community index can be 
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identified. This approach differs from discriminant analysis primarily in the nature of 
the dependent variable. In discriminant analysis, the dependent variable reflects cluster 
membership, based on overall species composition, while a community index can reflect 
more specific attributes of species composition in a quantitative manner. 

The community index used in this approach was a measure of dominance by a particular 
group of taxa (or guild). Guilds were defined by R-mode cluster analysis, which groups 
species based on similar spatial distributions (i.e., similar densities at each station). 
Densities were log transformed, as in the Q-mode cluster analysis of Section 2.2.2. An 
index of dominance by each guild was calculated at each station as follows: 

^log(Xi.l) 



S 

i=l 



where; S„ = number of taxa in the guild, 
o 

S = total number of taxa present, and 
X; = density of species i. 

The dominance index ranges from zero to one in magnitude. 

Some guilds consisted of recognized pollution tolerant taxa. Their index values were 
plotted on study area maps to illustrate association with known point sources. Bivariate 
correlations with sediment variables were determined, and stepwise multiple regression 
analysis was performed to find the best set of sediment variables for prediction of the 
community index. 

Some sediment variables were transformed prior to use m discriminant, correlation or 
multiple regression analysis, in order to improve their statistical distributions. Trace 
metals, pesticides, solvent extractables, total organic carbon (TOO, total Kjeldahl 
nitrogen (TKN) and phosphorus were log-transformed. An angular transformation was 
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1 ... 

used for particle size variables , and a square root transformation for loss on ignition 
(LOI). The pH and Eh distributions were not transformed. Water depth was included as a 
predictor variable, with log transformation. 

Water depth and Eh values for each station are listed in Table A5.1 (Appendix 5). Values 
for other sediment quality variables are available from the Water Resources Branch, 
Great Lakes Section. Mean values for all stations and parameters are provided in Table 
A5.2 (Appendix 5). A few stations missing data for pH, Eh, LOI and pesticides were 
characterized by study area mean values of these parameters prior to transformation. 



The particle size class variables used were: 

Class 1: > 1000 um diameter 

Class 2: ^5 to 1,000 um 

Class 3: < 'fJ um 

Silt 

Clay 

Silt + clay 
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3J RESULTS AND DISCUSSION 
M 1985 Survey 

Species and densities of benthic organisms found at stations in the St. Marys River 
(densities/m^) in 1985 are summarized in Table A4.1A (Appendix h). Raw data (species 
and numbers per sample) are provided in Table A4,1B (Appendix a). Physical, chemical 
and biological observations recorded at each station are presented in Table A5.1 
(Appendix 5). Species assemblages are described and compared with those found in 
earlier surveys in Section 3A. 

Detailed sediment quality results are available from the Great Lakes Section, \^ater 
Quality Branch, and will be reported in a forthcoming document. 

3.2 Indicator Species in the St. Marys River 

In the following subsections, reference is frequently made to the relative tolerance or 
intolerance of benthic species to polluted conditions. Clean-water indicators are 
reported as more prevalent at locations well removed from effluent discharges, while 
tolerant forms prevail in the apparently impacted areas found downstream of the major 
point sources. The pollution tolerance of benthic invertebrates, as outlined below, 
typically refers to tolerance of pollution by organic matter and associated iow dissolved 
oxygen conditions, rather then to tolerance of toxic metals and organic compounds. 

Mayflies such as Hexagenia , Ephemera and Caenis are generally considered to be 
intolerant (Roback, 197**); Hexagenia has been reported as rare or absent in sediments 
contaminated by oily substances in the St. Marys River (Hiltunen and Schloesser, 19S3). 
The tumbriculid Styiodrilus heringianus may also be considered a clean water indicator. 
This species prefers sandy to gravelly substrates such as found in much of the study area, 
but is typically absent or reduced in numbers in disturbed areas such as near urban 
centres (Nalepa and Thomas, 1976). Other species such as the tubificids Aulodnlus spp., 
Potamothrix moldaviensis and Spirosperma ferox are mesotrophic indicators, and 
Rhyacodrilus spp. is usually found in oligotrophic waters (Cook and Johnson, 197^*); thus, 
these species are probably intolerant of heavy organic pollution. Nematodes apparently 
reach their highest densities in mesotrophic regions, indicating a low tolerance for highly 
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polluted habitats, and their absence in highly polluted areas such as Toronto Harbour 
(Golini, 1979) may be related to their inability to tolerate low oxygen conditions for 
extended periods of time. The polychaete Manayunkia speclosa is also considered 
indicative of moderate organic pollution (Poe and Stefan, 1975), but intolerant of severe 
organic pollution (Mackie and Qadrie, 1971). Other mesotrophic indicators include the 
chironomids Polypedilum , Nanocladius , Psectrocladius , Dicrotendipes and 
Microtendipes . Gjll-breathing gastropods such as Valvata and Amnicola are common in 
the relatively unpolluted, shallow waters of outer Nipigon Bay (Lake Superior), but not in 
the inner bay where heavy organic pollution from pulp and paper wastewater is 
characteristic (Freitag et aL, 1973); therefore, these species also appear to be relatively 
intolerant. 

Among the pollution-tolerant forms, tubificids and chironomids are most note*'brthy. 
High densities of the tubificids Limnodrilus spp. (including immatures without capiUiform 
chaetae), Tubifex tubifex , Quistadrilus multisetosus and llyodrllus templetoni (both 
including immatures with capiUiform chaetae) are characteristic of areas showing 
organic enrichment throughout the Great Lakes (Cook and Johnson, 1974; Brinkhurst and 
Cook, 197'*). Similarly, the chironomids Procladius , Chironomus and Cricotopus are 
common in polluted conditions in the Great Lakes (Cook and Johnson, 1974). 

There is relatively little information available in the literature regarding the impact of 
heavy-metal stress on freshwater benthic communities or individual species. Some data 
has demonstrated that the oligochaete L^ hoffmeisten is very tolerant of high 
concentrations of some heavy metals (\l inner _et^. 1980) and is often the only species 
present (Wenste! et ^., 1977a) under such conditions. Nematodes are also indicated as 
tolerant of heavy metals, with some species occurring in high densities in the sand and 
gravel filter beds of sewage treatment plants (Pennak, 1978), where high metal 
accumulations would be expected. Similar data on the effects of persistent organic 
compounds on benthic organisms appear to be either extremely limited or unavailable. 

3»3 Cluster Analysis of Spatial Patterns in 1985 

Assignment of sampling stations to clusters, based on similarities in 1955 species 
composition, is presented in Table A 3.1 and Figure A3. 2 (Appendix 3). 
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Seven clusters were defined in 1985 (Figure 3.1) with three clusters (Nos. 2 to k) 
identified as including pollution-impacted stations. The stations in these three clusters 
were dominated by individual indicator species or pollution-tolerant groups of species, 
and were associated with known industrial outfalls or downstream depositional areas. 
The remaining four clusters (Nos. 1, 5, 6, and 7) generally described stations that were 
upstream of effluent discharges, near the U.S. shoreline or well downstream of 
discharges. Characterization of station clusters is summarized in Table 3.1. Biological 
characteristics presented in this table include species common to all stations within the 
cluster, mean number of taxa and mean total density. Physical and chemical features 
include substrate type, water depth, presence of macrophytes, presence of visible oil in 
the sediment and current speed. 

Clusters 2 to » 

The pollution-impacted clusters (2 to i*), include stations with benthic populations 
dominated by tubificid oligochaetes, along with nematodes and polychaetes. Most of the 
species present were at least moderately tolerant of the effects of organic enrichment. 
Cluster U included stations with considerable reductions in number of taxa and density, 
possibly as a result of toxicity from point-source discharge areas. These individual 
clusters are described in detail below. 

Cluster 2 

Cluster 2 consisted of five stations downstream of the eastern STP (SMD 5.0E-210, S.WD 
6.'tE-350, -f20, -1*90 and -570) and two other stations (Stations 113 at the north end of 
Lake Nicolet and 120 near the Algoma slip). These stations were characterized by a 
coarse sand substrate, with or without an overlying silt layer. Water depth was variable, 

« 

ranging from 2.5 to 1^ m, and current speed ranged from slight to moderate. 
Macrophytes were absent or very sparse at these locations, and oil was generally present, 
especially at SMD 6.*E-350. 

Characteristic taxa were nematodes, polychaetes, nemerteans, pea clams, immature 
tubificids with capiUiform chaetae (likely T. tubifex and/or L^ tempietom ), immature 
tubificids without capilliform chaetae (likely Limnodrilus spp.), S. ferox , and the 
lumbriculid worm S. heringianus . The mean number of taxa per station (23) was 
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Figure 3.1 

Relationships between 1985 Benthic 
Station Clusters Derived by 
Ward's Method 
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TABLE }.l: 



CHARACTERISTICS OF OnNTIIIC COMMUNITY CLUSTERS IN Tilt ST. MAUYS Kivru - 1')$^ 



Par.iineter 



CLUSTER NUMBER 
4 



COMMON TAXA 



Irnrndt. Tubilicids 
w/o chactae 

•Nematoda 

Pfocladius 

Be7zia 

' Manayunkia 

Spifospcrma 






•Immat. Tubilicids 

w/o cKaetae 

•Nematoda 

* Manayunkia 

Nemertea 

Stylodrilus 

Pisidium 

Spirospcrma 

Immature Tubificidi 

with chaetae 



•linmat. Tubihcids 
w/o chactae 
•Nematoda 
Nais variabilis 



•Ifnmat. Tubificids 
w/o chaetae 
•Nematoda 
•Chironomidae 



Immat. Tubificids 

w/o chaetae 

•Nematoda 

Procladiiis 

Cricotopus 

Bezzia 



Immat. Tubilicids 
w/o chaetae 
• Nematoda 
Martayunkia 
Pisidium 
Amnicota limosa 
Valvaia sincera 
Procladius 
Polypeditum 
•Chironomidac 



Immat. Tubificids 

w/o chaetae 

•Nematoda 

He U soma 

Procladius 

Chironomus 

Polypedilurn 

HyallelJ 

Hydracarina 

•Chironomidac 



WEAN NO.TAXA 27 

MEAN TOTAL 

DENSITY (No/m^) 128,000 



SUBSTRATE 



Drown Sill Over 
Sand/Clay 



WATER DEPTH (m) 2 13.7 

MACROPHYTES Variable 

VIMIU E OIL Absent-Slight 

CUIiRI NT Usually None 
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192,000 

Coarse Sand 
w or w/u Sill 



13 

2)9,000 
Organic Silt 



2.5-l<» 1-16 

Absent oi Sparic Variable 

Absent Very Strong Stitht-Veiy Strong 

Slight-Moderate None 



12 

7 1 .000 
Sill 



39 



56,000 



Variable 

Silty-("ojf se Sand 



1.5 8.5 1-5.5 

Usually Absent Variable 

Absent-Very Strong Absent-Slight 

None-Strong Noiu-Moderale 



40 

165,000 

Sill Over Sand 

J. 5-10. 5 
Present 
Absent-Slight 
None-Slight 



32 

201,000 

Brown Silt Over 
Sand/Clay 

JI2 

Usually Present 

A blent -Very Strong 

None-Slight 



noiiiiriant taxa 



2 

relatively low, while the mean total density was relatively high (192,000 organisms/nn ) 

in comparison with the apparently unimpacted group of station clusters (Nos. 1, 5, 6 and 
7). Most of the common species among Cluster 2 stations were indicative of mesotrophic 
conditions, with more sensitive groups such as gastropods, isopods, ephemeropterans and 
trichopterans either substantially reduced in density or absent. On the basis of the 
relatively low numbers of taxa and high total densities relative to control areas upstream 
of Sault Ste. Marie and along the U.S. shoreline, stations included in Cluster 2 were 
considered to be slightly impaired in terms of pollution status. 

Cluster 3 

Stations included in Cluster 3 were situated near the Canadian shore of the St. Marys 
River, downstream of the main industrial outfalls (SMD 1.2-8'^0, SMD 2^-950); 
downstream of the eastern STP (SMD 5.0E-370 and -6W); in Little Lake George (S6, 87, 
88); and in Lake George (100 and 102). All of these stations were depositional areas 
characterized by a silt substrate, no obvious current and the presence of an oily sheen on 
the sediment surface. This oily sheen ranged from slight to heavy, with the heaviest 
accumulations noted at two stations downstream of the industrial discharge - SMD 1.2- 
8^*0 and SMD 2.6-950 (based on interpretation in the field). Macrophytes (Potamogeton 
sp. and Elodea ) were present at these two stations and Stations 87 and 88 in Little Lake 
George. 

Characteristic organisms at Cluster 3 stations consisted of nematodes, immature 
tubificids without capiUiform chaetae (likely Limnodrilus spp.) and the naidid oligochaete 
Nais variabilis . Other frequently abundant taxa at these stations included the tubificids 
q. multisetosus , L^ templetoni and/or T. tubifex (both including immatures with 
capilliform chaetae), and chironomid species such as Chironomus sp., Dicrotendipes , 
Polypedilum , Cricotopus , Procladius and Ablabesmyia . The snail Helisoma sp. was also 
frequently present. All of these species (except Nais ) have been identified as tolerant of 
severe organic enrichment. In addition, the relative abundance of the tubiflcid fraction 
was substantially greater than observed at Cluster 2 stations. The presence of 
nematodes suggests, however, that the conditions were not as impaired as those in 
Toronto Harbour where nematodes were reported as nearly absent in 1977 (Golini, 
1979). The largest Nais variabilis populations were generally associated with areas 
supporting macrophyte growth. Of particular note was the low density or absence of the 

- 26 - 



polychaete M. speciosa at Cluster 3 stations, in contrast to its abundance at stations in 
the closely related Cluster 2 {Figure 3.1). The absence of M. speciosa suggests greater 
impairment of water and sediment quality than occurred at stations in Cluster 2. This 
greater degree of impairment is reflected in a much lower mean number of taxa (15) 
relative to Cluster 2 and a much greater mean total density (259,000 organtsms/m ). 
Both of these features are typical of organic enrichment. On this basis, as well as the 
species assemblage, Cluster 3 stations were considered moderately impaired. 

Cluster 4 

Cluster li included eight stations in the immediate vicinity of the Algoma Steel slip (SMU 
1.0-2130, SMU 0.5-300, Stations 115, 116, 117, 118, 119 and 127), two stations near the 
Canadian shore below the rapids (SMD O.S-600 and Station 107) and one station 
immediately downstream of the eastern STP {SMD 5.0E-520). Physical environmental 
conditions described at stations within this cluster group were variable. Substrate 
conditions varied from silt (SMD 5.0E-520) to hard pack clay (Station 116), with water 
depth ranging from 1.5 to 8,5 m. There were generally no macrophytes, and the presence 
of an oily sheen varied from absent to heavy. Current conditions were also variable, with 
descriptions ranging from no obvious current to strong current. 

The cluster was defined biologically by the presence of immature tubificids wirhout 
capilUform chaetae (likely Limnodrilus spp.). Nematodes were also present at a majority 
of stations. Other common species included the tubificid Q. multisetosus and the 
chironomids Procladius and Cryptochironomus . Few other taxa were common at these 
stations. All of the species described above are tolerant of organic pollution. The mean 
number of taxa per station in Cluster i* was only 12 (the lowest of all clusters), with a 
mean total density of 71,000 organisms/m^. The variation in total densities at stations 
within Cluster U is very large, with very high densities occurring at organically enriched 
areas, and very low densities at some locations where effluent discharges may be 
exerting a toxic effect. For example, the mean total density of organisms at Station 1 16 
in the Algoma slip was only 88/m , suggesting a toxic inhibition of benthic community 
development, perhaps in combination with poor habitat conditions (the substrate was a 
rather hard, uniform clay with very little overlying silt). At Station SMD 5.0E-520 near 
the STP, severe organic enrichment is implied by the extremely high mean total density 
of 345,000/m^. One feature common to all stations, however, is a low number of taxa 
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relative to stations in ai! other dusters, and a general predominance of pollution-tolerant 
forms. 

Biologically unusual stations within Cluster t included Stations 127 and SMD 0.8-600. At 
Station 127 in the presence of a silt substrate, nematode densities reached 220,000 
organisms/m,^ or 70% of the population. The remainder consisted primarily of 
L. hoffmeisteri (including immatures) and ^ multisetosus . At Station SMD 0.8-600, 
coelenterates (Hydra sp.) accounted for 51% of the population, with the nemertean 
Prostoma rubrum and the oligochaete family Enchytraeidae accounting for the majority 
of the remainder. Both of these stations are near effluent discharge sources, and the 
biological assemblages probably reflect the Influence of these sources. Unfortunately, 
information on the pollution tolerances or indicator value of most of these organisms is 
scant. 

Clusters 1, 5, 6, and 7 

The remaining four clusters (1, 5, 6 and 7) included stations which were either upstream 
of the discharges, near the U.S. shore or too far downstream to be substantially 
affected. As such, they were characterized by benthic associations which were low in 
tubificid oligochaetes and are considered to be relatively unimpaired. 

Cluster 1 

Cluster 1 included the majority (14 of 22) of stations located in Lake George. Other 
stations in this cluster included one at SMD 5.0E-60 and Station 85 in Little Lake 
George. All of these stations were characterized by a brown silt substrate with either 
sand or clay underneath. Depth was variable, ranging from 2 to 14 m, and current was 
generally non-existent, with three exceptions where there was a slight current present. 
The presence of macrophytes was variable, while oil tended to be visibly absent or 
present only in trace amounts. 

Stations in Cluster 1 were characterized by five primary taxa, which included nematodes, 
polychaetes (M. speciosa ), immature tubificids without capiiliform chaetae (likely 
Limnodnlus spp.), the chironomid Procladius and the ceratopogonid Bezzia . Other 
species common to the majority of the stations in this cluster were the pea clam 
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Pisidium spp., the tubificid 5. ferox , immature tubificids with capilliform chaetae (likely 
T. tubifex and/or I. templetoni ) and the chironomids Cryptochironomus and 
Cryptotendipes . 

The mean number of taxa per station in Cluster 1 was 27 (range 12 to 39), with a mean 
Total density of 128,000 organisms/m^. High organism densities and the presence of non- 
toierant organisms at most of the stations within the cluster, suggest that the region 
including Cluster 1 stations was generally unimpaired overall, with some degree of 
organic enrichment. Variation within the cluster does exist, however, with stations at 
the north end of Lake George (Stations 96, 97, 98, iCf, 105 and 106) tending toward 
moderate impairment, displaying unusually high densities of either nematodes or 
polychaetes and having an oily sheen on the surface of the substrate. This area likely 
represents a depositional zone from the St. Marys River inflow. 

Cluster 5 

Stations included in Cluster 5 were situated downstream of Whitefish Bay (Station 126^ 
2.5 km above the St. Marys falls (Stations SMU 2.5-1830, -1300 and 123, 12^*) on the 
Canadian shore upstream of the Algoma slip (Station 121), opposite the Algoma slip (SMU 
0.5-200), and on the northeast side of Lake George (Station 99). Station characteristics 
were variable for all recorded parameters. Substrates ranged from siity-sand to coarse 
sand with a relatively shallow water depth range of 1 to 5.5 m. Visible oil was generally 
absent or at the most slight, and macrophytes were present or absent. Perceived current 
was absent in most cases, with only Station 121 and SMD 0.5-200 having a moderate 
current. 

Characteristic taxa were Nematoda, immature Tubificidae (without capilliform chaetae), 
the ceratopogonid Bezzia and the chironomids Procladius and Cricotopus . Other 
frequently occurring taxa included non-tolerant forms such as the ephemeropterans 
Hexagenia , Ephemera , Paraleptophlebia , Caenls and Stenonema ; trichopterans such as 
Molanna and the snail Valvata . Tubificids tended to be a minor portion of this 
community, with chironomids forming the majority of the assemblage. While most of 
these tended to be recognized as facultative forms (i.e., Procladius and Cricotopus ), the 
high mean number of taxa (39), the low mean total density (56,000/m ), and the location 
of these stations either upstream or well downstream of known point sources in the river 
suggest that this population is representative of a relatively unimpaired community. 
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Cluster 6 

Cluster 6 was the second largest cluster (after Cluster 1), including It* stations. These 
stations were located in Leigh Bay (SMU 2.5-2k7Q), near the U.S. shore below the St. 
Marys Falls (SMD 1.2-15 and SMD 2.6-W), in Lake Nicolette (Stations 108, 109, 110, 111, 
112 and l\k) and in Lake George (Stations 83, 90, 91, 9^* and 103). Characteristics of 
these locations were a substrate community of silt over sand, a water depth ranging from 
1.5 to 10.5 m, a local macrophyte community (Chara , Potamegeton , Elodea and Isoetes ), 
an absence or very slight trace of oil on the substrate surface, and generally no current, 
indicating potential depositional areas. Essentially all locations were either upstream or 
across the river from point sources, or were well removed downstream (Lake George). 

Characteristic taxa included the ubiquitous nematodes and immature tubificids (vJ/lthout 
capilliform chaetae), the chironomids Procladius and Polypedilum , the polychaete M. 
speciosa , pea clams, and the snails V. smcera and Amnicola limosa . Other commonly 
occurring taxa included the isopods Asellus and Lirceus , the amphipod Hyalella azteca 
and the ephemeropterans Hexagenia and Ephemera . This assemblage is a mixture of 
facultative and non-tolerant organisms. The mean number of taxa at each station in this 
cluster was the highest of all seven at kQ, with a mean total density of 163,000 
organisms/m^. As in Cluster 5, chironomid taxa dominated the benthic populations at 
the stations included in Cluster 6. This cluster is represented in situations that may be 
classified as relatively unimpaired, based on the high number of taxa, the presence of 
intolerant forms, and remote location of these stations from point sources of pollution. 

Cluster 7 

Cluster 7 is the smallest of all the clusters, and included only five stations. These were 
located upstream of the St. Marys Falls (Station 122, 125 and SMU 1.8-1690), downstreair. 
of the eastern STP (SMD 6.^E-2S0), and in Little Lake George (Station 89). Stations m 
this cluster were characterized by a brown silt over sand/clay substrate, water depths of 
3 to 12 m, a general presence of macrophytes and very little if any perceptible current. 
Visible oil on the substrate surface was variable, and ranged from none to very abundant 
(SMD 6.fE-280). 
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Taxa common to all stations included nematodes, immature tubificids (without 
capilliform chaetae), the chironomids Chironomus , Procladius and Polypedilum . 
Hydracarina (water mites), the amphipod hL azteca and the snail Helisoma . Other 
common taxa included the isopods Asellus and Lirceus and the chironomids 
Psectrocladius , Polypedilum , Ablabesmyia and Procladius . As noted before, most of 
these taxa are facultative or tolerant forms. However, the mean number of taxa within 
this cluster remains relatively high at 32, with a mean total density of 201,000 
organisms/m , the majority of these being chironomids. Tubificids were found in 
relatively low densities. In general, stations within this cluster may be classified as 
unimpaired by virtue of the relatively high mean number of taxa. Exceptions include 
some stations identified as having a high oil content in the sediments, although these 
locations supported benthic faunas that were generally similar to those having lower oil 
levels within the same cluster. 

3mA Tempered Comparisons of Spatial Patterns 

Results of the 1985 cluster analysis are summarized in Figure 3.2A and 3.2B in terms of 
the following pollution impairment zones: 

1. severe: 

a) extreme tubificid dominance (i.e. , L. hoffmeisteri and immatures 
without capilliform chaetae), tolerant chironomids, low numbers of 
taxa and high total densities, or 

b) communities with either very low total densities and low numbers of 
taxa, and/or high densities of nematodes with few other taxa (Station 
Cluster k); 

2. moderate: tubificid dominance with high densities of nematodes and 
facultative chironomids, absence of polychaete worms, reduced numbers of 
taxa and high total densities (Station Cluster 3); 

3. slight: nematode and polychaete dominance with moderate densities of 
tubificids and some non-tolerant groups present (Station Cluster 2); and 

li. unimpaired: communties tending towards chironomid dominance, with 
several non-tolerant groups (e.g., ephemeropterans and trichopterans) 
present, low tubificid densities and high numbers of taxa (Station Clusters 
1,5,6 and 7). 
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Figure 3.2a 
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Figure 3.2b 



Distribution and Zones 
of Impairment of 
Benthic Fauna-1985 



The spatial distribution of these four zones, based on assignment of each duster of 
stations to one of the zones, is used to infer any major changes in benthic community 
status between years. 

In comparisons among surveys, methodological differences must be considered. As noted 
in Section 2, increased sorting efficiency (due in part to smaller mesh size and use of a 
microscope) contributed to higher total densities in 1985 than in earlier surveys. Total 
organism densities in 1985 ranged from 88 to 591,000 organisms/m , compared to a range 
of ICf to ifO.OOO/m^ in 1983, a difference of up to one order of magnitude. 

In addition, principal components analysis (Section 2 and Appendix 2) suggested numerous 
faunal shifts between 1983 and 1985. Two particularly notable shifts involved nematodes 
and polychaetes which were rare or absent in all previous sureys, and abundant and often 
dominant in 1985. These organisms are small enough to be influenced by mesh size or 
sorting efficiency (Golini, 1979). Other small organisms, which occurred in higher 
densities in 1985, included naidid oligochaetes and early instar chironomids. However, 
other shifts such as the increased densities of Hexagenia , a comparatively large mayfly 
nymph, in 1985, were probably not related to mesh size or sorting technique. 

Total organism densities reported at some locations in 1985 are very high ( >100,000/m ) 
and comparable to the maximum densities reported in the most heavily polluted harbours 
of the Great Lakes (Cook and Johnson, 197ft). These high densities can, to some extent, 
be accounted for by the large numbers of very small organisms retained by a 200 um 
mesh. In contrast, most benthic surveys of the Great Lakes have been carried out using a 
U.S. #30 (about 500 um) mesh, which would not retain as many organisms. Greater 
research into the environmental tolerances of many of these small but abundant taxa 
(nematodes, naidids) would facilitate interpretation of their significance in benthic 
communities and pollution surveys. 

In spite of the differences between surveys, inferred zones of impairment in 1985 were 
similar to those reported in 1983 (Figure A 1.3). Both surveys indicated a zone of severe 
impairment below the St. Marys Falls near the Canadian shore at SMD 0.8. In 1985, 
however, this zone did not include SMD l.2-8ft0, suggesting some improvement here since 
1983. This change was primarily due to the presence in large numbers of the nemertean, 
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Prostoma rubrum , which is typically found in weli-oxygenated, littoral standing water 
and cannot tolerate low oxygen conditions (Pennak, 1978). 

The 1985 survey indicated severe impairnnent, possibly resulting from toxicity, in the 
vicinity of the Algoma slip. In 1983, this area was identified as moderately impaired, 
based on 3 stations located outside the slip. The greater nunrvber of stations in closer 
proximity to the slip probably accounts for the change in status of this area in 1985 
relative to 1983. Species assemblages at the stations common to both surveys had not 
changed markedly. 

Both surveys also suggested a zone of moderate impairment which extends downstream 
from St. Marys Falls, close to the Canadian shore, to the vicinity of 5MD 2.6. Another 
zone of moderate impairment was associated with the easterly Sault Ste. Marifr STP, 
extending to within 370m of the U.S. shore. Differences from the 1983 pattern included 
the identification of SMD 5.0E-520 and station 107 as severely impaired in 1985. 
Analysis of the detailed species lists from both years substantiates the differentiation of 
SMD 5.0E-520 from adjacent stations SMD 5.OE-640 and 370 on the basis of a more 
complete tubificid dominance (primarily L^ hoffmeisteri ) and only 3 other taxa. The 
differentiation of station 107 from nearby stations was related to a lower total density 
and fewer taxa. 

In 1985, a zone of slight impairment was evident within the Lake George Channel 
discharging into Lake George. In 1985, this zone extended from SMD 5.0E-2I0 
downstream past transect SMD S.'fE (except SMD 6.<^E-210 near the U.S. shore), into 
Little Lake George. In the early previous surveys (1968 and 1973), this section was 
designated as moderately impacted. The apparent improvement here is based on the 
presence of taxa such as S^ heringianus , P. rubrum as well as on a slight increase in taxa 
and reduction in densities relative to the moderate zones. The absence of an impairment 
designation in this portion of Little Lake George in 1983 is attributed to the lack of 
sampling in this area during the 1983 survey. Another station described as slightly 
impacted in 1985 includes Station 120 just outside the severe impact zone near ;he 
Algoma slip. 

A zone of moderate impairment at Stations 86, 87 and 88 was identified in 1985 but not 
1983 when these locations were not sampled. Species assemblages were similar here to 
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those reported in 1973 and were characterized by high tubificid densities (primarily L. 
hoffmeisteri and low nunnbers of taxa). The extent of this zone is sirr.ilar to that 
indicated in 1973 (Hamdy et aL, 1978). 

A zone of moderate impairment was also identified in Lake George, at the two deep 
water stations. The sampling in Lake George in 1985 was the first time the lake has been 
extensively sampled, apart from a study in 1975 by Hiltunen and Schioesser (1983). 

Hiltunen and Schioesser (1983) related the distribution of the ephemeropteran Hexagenia 
to the presence or absence of visible oily residues in sediments in the St. Marys system 
downstream of the STP . Substrates collected from Lake George in 1985 were 
contaminated by oily substances at most stations between the STP and upper Lake 
George, similar to the pattern indicated by Hiltunen and Schioesser. The 1985 results 
suggest that the extent of sediment contamination by oU in the Lake George Channel has 
not changed from that observed in 1975 (Figure 3.3). Oil was absent in the Lake Nicolet 
Channel in 1985, and was very limited in distribution in the more extensive study of this 
area in 1975 (Figure 3.3). The distribution of Hexagenia in 1985 has remained similar to 
that observed in 1975. In 1985, Hexagenia was absent from sediments in 27 of 32 stations 
analyzed in the St. Marys River, Lake Nicolet, Lake George Channel, Little Lake George 
and Lake George where even slight amounts of visible oil were present. In the remaining 
19 stations where visible oil was absent, 16 had Hexagenia populations present (Figure 
3.3). At two of the three remaining stations where Hexagenia was absent, the bottom 
substrate was primarily sand which deters Hexagenia inhabitation (Edmunds, M^ 1976). 
Where similar stations were sampled, densities of Hexagenia were usually higher in 1985 
than in 1975 (Figure 3.3). Whether this difference represents a real temporal trend or is 
due to differences in sampling is not known. 



This does not imply that the STP is the source of the oil. Oil also occurred in 
abundance at stations upstream of the STP in 1985. 
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FIGURE 3.3 DISTRIBUTION ON HEXAGENIA NYMPHS AND 

VISIBLE OIL IN THE ST. MARYS RIVER SEDIMENTS 
IN 1975 AND 1985 
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3.5 Sediment Quality Relationships of Benthic Communities 

3.5.1 Distributions of Sediment Quality Variables 

Prior to transformation, 29 of 33 sediment parameters were significantly skewed, 27 
showed significant (p < 0.05) kurtosis, and none were normal in distribution. After 
transformation, 13 sediment parameters were still skewed, 15 showed kurtosis, and 18 
were normal in distribution. In all cases, the selected transformation improved the 
distribution, by reducing the absolute magnitude of skewness and kurtosis measures 
(Appendix 5, Table A5.2). 

It was not possible to normalize any of the pesticide parameters due to the large number 
of observations below the detection limit (these were all considered equal- to the 
detection limit for statistical analysis). Similarly, pH was unlikely to be normalized due 
to the large number of missing observations replaced by the mean. Subsequent analyses 
were performed with and without the non-normal parameters; however, results based on 
the full data set were most interpretable, and these are presented in the following 
sections. 

3.5.2 Discriminant Analysis 

Figure 3.4 illustrates separation in discriminant space of the seven station groups defined 
by cluster analysis in Section 3.3. Stations in the same cluster (with the same numeric 
label in Figure 3.4) have similar overall species compositions. Clusters 2, 3 and 4 were 
considered to represent impacted groups, based on the recognized pollution tolerance of 
their characteristic species (Table 3.1, Section 3.3). 

Station Clusters 3 and 4 are distinguished by high sediment concentrations of iron and 
zinc in relation to particle size (DFl in Figure 3.4). Cluster 2 sediments are 
distinguished on the same discriminant function by a large fraction of particulate 
material (45 to 1,000 um), but lower concentrations of these metals. Relationships 
between bulk sediment chemistry as examined here, and benthic community response, are 
poorly understood, and may be strongly influenced by the fractions of sediment 
contaminants that are loosely associated with the sediment and therefore are more 
bioavailable. Nonetheless, geometric mean iron concentrations of 3.1% and 6.3% and 
zinc concentrations of 174 and 257 ug/g at stations defined by Clusters 3 and 4, 
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FIGURE 3.4 

Station Clusters Plotted on First Two 
Discriminant Functions 
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respectively, are considerably higher than the MOE dredge disposal guidelines of 1% for 
iron and 100 ug/g for zinc (Persaud and Wilkens, 1976), indicating that these key 
elements in DFJ are relatively high in concentration at these impacted locations; thus 
their influence on cluster separation appears plausible. Cluster 3 sediments are 
distinguished from those of Cluster i* by higher concentrations of PCB and less 
particulate material of size class 2 (DF2), in addition to their metal content. The first 
three discriminant functions were statistically significant (p < 0.05), with the first two 
collectively accounting for 72% of variance in cluster membership (Table 3.2). 

Station Clusters 1, 5, 6 and 7 were not readily distinguishable from each other on the 
basis of the first two discriminant functions, but were reasonably well distinguished as a 
group from each of the three impacted station clusters. With the other functions, 
primarily DF3 (a particle size gradient), some discrimination within this group was 
achieved. Cluster 1 had the most fine material { < <t5 urn), while Cluster 3 had the least. 
Three stations within this four-cluster group were misclassified by the full set of 
discriminant functions (i.e., predicted on the basis of sediment characteristics to be ir. 
the wrong cluster). The overall percentage of stations correctly classified was 95.7% 
(Table 3.3). 

Standardized discriminant function coefficients are tabulated in Appendix 6, Table 
A6.1. Sediment variables with coefficients greater than 1.9 in absolute value are listed 
on the DPI and DF2 axes of Figure 3A. Means for these variables in each station cluster 
(after transformation) are shown in Table 3.'*. Discriminant function scores for each 
station are tabulated in Appendix 6, Table A6.2. 

It was not possible to perform Box's test of multivariate homogeneity of variance with 
this analysis, because there were too few stations in each cluster in relation to the 
number of sediment variables. This results in non-singular group covariance matrices 
which the test cannot accommodate. 

Varimax rotation of the discriminant functions was attempted (Kaiser, 195S). However. 
the unrotated solution was more interpretabie. Green (1979) notes that this is often the 
case in ecological applications. The rotated solution is not presented. 
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TABLE 3.2: SlGiNIFICANCE OF DISCRIMINANT FUNCTIONS OF SEDIMENT 

CHARACTERISTICS 



Eigen- Percent of Cumulative SJRnificance After Removal 

Function value Variance Percent Chi-squared dl p 



10.it8398 1*2.10 03.10 278,0{f 170 Q.OOQO 

7.12867 29.30 72M 177M 132 0.0051 

2.87730 11.83 8i>.23 112,i*2 96 0.1209 

1.90372 7.83 92.06 61.25 62 0.5031 

l.if9856 6.16 98.22 17.23 30 0.9689 

0.43375 1.78 100.00 - - - 



Standardized discrinninant function coefficients for each sediment variable, and 
discriminant scores on each function for each station, are listed in Appendix 6. 
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TABLE 3.3: PREDICTION OF STATION CLUSTER MEMBERSHIP FROM 

DISCRIMINANT FUNCTIONS OF SEDIMENT CHARACTERISTICS 



Actual No. of 

Cluster Stations 



Predicted Cluster Membership 



1 


16 


15 


^ 


2 


7 


- 


7 


3 


f 


- 


- 


4 


11 


- 


- 


5 


S 


- 


•w 


6 


14 


- 


- 


7 


5 







11 







Percent 


6 


7 


Correct 


1 


- 


93.8 


- 


- 


lOO.O 


• 


• 


100.0 


- 


- 


100.0 


- 


- 


100.0 


13 


1 


92.9 



so.o 



TOTAL 



70 



95.7 
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TABLE 3.4: STATION CLUSTER MEANS OF SELECTED^ SEDIMENT 

CHARACTERISTICS 



Sediment 




Station 


Cluster M 


ean (after transform 


ation)^ 




Variable^ Cluster: 


1 


2 


3 


k 


5 


6 


7 


Iron 


*.38 


3.90 


kAS 


4.80 


3.96 


3.91 


4.35 


Zinc 


2.09 


lAi 


2.214 


2.41 


1.35 


1.47 


1.99 


Size 2 ik5-i,QQ0 urn) 


0.66 


1.23 


0.71 


0.93 


1.15 


1.C7 


0.76 


Size 3 i<^5 um) 


0.91 


0.30 


0.85 


0.60 


0.36 


0.49 


0.80 


Aluminum 


4.06 


3M 


3.96 


3. 89 


3.57 


3.64 


3.94 


PCBs 


1.33 


1.30 


1.57 


1.40 


1.30 


1.38 


1.36 


Endosulphan II 


0.61 


0.60 


0.65 


0.77 


0.60 


0.6: 


0.60 


Silt 


0.92 


0.32 


0.87 


0.52 


0.43 


0.57 


0.S5 



' Sediment variables shown are the main contributors to discriminant functions I and 2, 
with standardized DF coefficients exceeding 1.9 in absolute value. 

See Table A5.2 for original units and transformation. 
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A stepwise discriminant analysis was also performed, using four groups of stations 
{Clusters 1, 5, 6 and 7 combined). In this analysis, a sediment variable is removed from 
the predictor set (if it ceases to make a significant contribution), or added to the 
predictor set (if it significantly increases discrimination between groups) at each step. 
Green (1979) recommends against the stepwise approach, since the discriminant functions 
rather than the original variables are assumed to be controlling species distribution; 
however, it is one way of reducing the number of original variables to the point that 
Box's test can be performed. The stepwise solution is illustrated in Appendix 6, Figure 
A6.1. 

3.5.3 Community Index Analysis 

Benthic species clusters, or guilds, based on R-mode cluster analysis, are illustrated in 
Figure 3.5. Similarities between guilds with respect to their spatial distributions are 
indicated by the rescaled distance values at which they combine. Distance is an inverse 
measure of similarity. Five guilds are evident at an interpretation level of 13.5 on the 
distance scale. These are labelled A to E on the dendrogram. 

Guilds B and C are comprised primarily of pollution tolerant taxa (Figure 3.5). Guild B 
consists of Nematoda and immature Tubificidae without hair setae. The immatures are 
most likely L. hoffmeisteri based on the relative abundance of adults. Station clusters 2, 
3 and k were all characterized by dominance of Guild B taxa, and were identified as 
pollution-impacted (Section 3.3). Guild C includes both pollution-tolerant species, such 
as Limnodrilus spp., Q. multisetosus and Chironomus , and some species that are less 
tolerant of organic pollution, such as Aulodrilus spp., suggesting that Guild C occurs in 
areas that are less impaired than areas inhabited by Guild B. 

Species Guilds D and E are described by some relatively intolerant and mesotrophic taxa, 
such as Hexagenia (Guild D) and Polypedilum (Guild E), with relatively few eutrophic 
indicators. Guild D is probably more mdicative of unimpacted conditions than is Guild E, 
based on the occurrence of the intolerant HexaRenia in the former, and greater 
representation of forms tolerant of organic enrichment in the latter (e.g., Tubificidae, 
Procladius , Manayunkia ). High index values for these guilds, particularly for Guild D, 
tend to occur at stations assigned to Station Clusters 1, 5, 6 and 7, which are described 
as relatively unimpaired (Section 3.3). 
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FIGURE 3.5 Species Clusters (Guilds) Based on 
Concordance of Spatial Distributions 
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Species Guild A includes a very large number of species ( >100), most of which tend to be 
uncommon at most locations. Because most of the dominant forms and many of the key 
Indicator species occurred in Guilds B to E, interpretation of the pollution indicator value 
of Guild A was not attempted. 

The community indexes representing dominance of each Guild B, C, D and E were 
examined for correlation with sediment characteristics. Indexes for each station are 
listed in Appendix 6, Table A6.3. Coefficients of correlation between dominance indexes 
and sediment parameters are listed in Table 3.5. 

Guild B and Guild C dominance indexes were each positively correlated with various 
heavy metals, solvent extractables and TOG (p < 0.01). Guild B dominance was also 
positively correlated with several pesticides, while Guild C dominance was also 
correlated with fine particulates, phosphorus, LOI and TKN. Guild B dominance was 
negatively correlated with Eh, indicating association with reducing environments. These 
correlations support the contention that dominance by Guilds B and C is representative 
of impaired conditions. 

Guild E dominance was positively correlated only with depth (p < 0.0 1), and negatively 
correlated only with ppDDE, suggesting that these taxa are not strongly indicative of 
pollution status. Guild D dominance was negatively associated with various heavy 
metals, endrin, solvent extractables, phosphorus, LOI, TOG and TKN. Thus, this group 
represents a relatively clean-water faunal assemblage, both on the basis of biological 
indicators and negative correlations with sediment-associated contaminants and organic 
matter. 

Based on correlations with sediment characteristics, Guild B dominance appears to be the 
best general indicator of toxic contaminants (both metals and pesticides). This guild is 
also defined by two of the most abundant indicator taxa (Nematoda and Limnodnlus 
spp.), thus facilitating further interpretation of spatial trends in environmental quality. 
This index was plotted on a map of the study area to illustrate associations with known 
point sources. Figure 3.6 shows the highest index values in the vicinity of the Algoma 
slip and the Sault Ste. Marie, Ontario STP. Values greater than 0,25 in these areas were 
considered indicative of severe impact. Zones of moderate impact (index 0.18 to 0.25) 
were also apparent in a depositionai area downstream of St. Marys Falls on the Canadian 
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TABLE 3.5: CORRELATIONS OF GUILD DOMINANCE INDEXES WITH SEDIMENT 

CHARACTERISTICS 



Correlation Coefficient 



Sediment 
Variable 



Guild: 



Water Depth 

Aluminum 

Arsenic 

Cadmium 

Cobalt 

Chromium 

Copper 

Iron 

Mercury 

Manganese 

Nickel 

Lead 

Selenium 

Zinc 

Silt Fraction 

Clay Fraction 

Silt + Clay 

Coarse (> 1,000 urn) 

Medium ('*5- 1,000 urn) 

Fine (< 45 um) 

Dieldrin 

Methoxychlor 

Eh 

Endrin 

Endosulphan II 

Endosulphan I 

Heptachlor Epoxide 

PCS 

ppDDE 

Loss on Ignition (LOI) 

Solvent Extractables 

Total Organic Carbon (TOO 

Total Kjeldahl Nitrogen (TKN) 

pH 

Total Phosphorus 



-0.0353 


0.1052 


0.3572* 


-0.3148* 


O.'HOS* 


0.1489 


0.1470 


-0.3976* 


0.1786 


0.4349* 


0.0442 


-0.3459* 


0.'*153* 


0.0164 


0.1782 


-0.2315 


0.21*13 


-0.0832 


0.0729 


-0.1733 


0.4104* 


0.2491 


-0.0816 


-0.3543* 


0.4740* 


0.3016* 


-0.0059 


-0.4349* 


0.2189 


0.4115* 


0.0233 


-0.'#154* 


0.5089* 


0.3946* 


0.1442 


-0.5114* 


0.0916 


0.3284* 


0.0507 


-0.2278 


0.3443* 


0.2517 


-0.0334 


-0.3348* 


0.3660* 


0.3763* 


0.0863 


-0.4392* 


0.5006* 


0.2153 


0.0430 


-0.3899* 


0.3641* 


0.4413* 


0.1321 


-0.4800* 


0.4941* 


-0.0212 


0.0963 


-0.1950 


0.2536 


0.1128 


0.1272 


-0.3189* 


0.4627* 


0.0075 


0.1048 


-0.2629 


-0.2372 


0.1407 


-0.1760 


0.0951 


-0.4136* 


-0.0625 


-0.1416 


0.3108* 


0.4324* 


0.0594 


0.1459 


-0.3082* 


0.1671 


0.0492 


-0.2389 


0.02S9 


0.0786 


0.5862* 


-0.1085 


-0.3805* 


-0.2797 


-0.3907* 


-0.0001 


0.2913 


-0.2382 


0.4975* 


-0.2562 


-0.0060 


0.0146 


0.7105* 


-0.1828 


-0.3171* 


0.1488 


0.5714* 


-0.1884 


-0.3189* 


-0.0628 


0.6914* 


-0.2753 


-0.2024 


0.2820 


0.1240 


-0.2842 


-0.1055 


0.3820* 


0.1367 


-0.3347* 


-0.1156 


0.4793* 


0.2855 


0.0098 


-C.3S15* 


0.5910* 


0.32S4* 


-0.0380 


-0.4236* 


0.3511* 


0.5807* 


-0.0831 


-0.4336* 


0.5408* 


0.0381 


0.0193 


-0.1888 


0.2294 


0.0664 


-0.2340 


0.1266 


0.5311* 


0.0668 


0.1068 


-0.3562* 



» Indicates coefficients differ significantly from zero (p< 0.0 i). 
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shore, downstream of the Sault Ste. Marie, Ontario STP to Little Lake George, and in a 
deep water depositional area of Lake George near the inlet. Mean sediment 
characteristics associated with each zone are listed in Table 3.6. 

Impact zones defined on the basis of Guild B dominance (Figure 3.6) closely resemble the 
zones previously defined in terms of station clusters 2, 3 and if (Figure 3.2, Section 3.3). 
However, they suggest a slightly greater degree and extent of impact in Little Lake 
George and Lake George than discerned by station cluster distributions. 

The best multivariate equation for prediction of Guild B dominance from sediment 
characteristics includes significant (p < 0.05) positive contributions from mercury, 
heptachlor epoxide, clay and TOG (Table 3.7). High concentrations of these constituents 
tend to increase the predicted dominance index. The equation includes significant 
negative contributions from arsenic, chromium and LOL The complete equation explains 
72% of variation in the dominance index (adjusted R = 0.72). 

The sediment characteristics identified by the multiple regression analysis as the best 
quantitative predictors of Guild B dominance differ from those identified by discriminant 
analysis as the key contributors to qualitative discrimination between station clusters 
defined on the basis of overall species composition. This is not surprising in view of the 
different dependent variables involved (indicator species dominance vs. station cluster 
membership). Nevertheless, the concordance of impact zones derived from the 
dominance index and station cluster patterns strengthens the finding of continued 
existence of tolerant benthic communities downstream of the Algoma Steel, St. Marys 
Paper and Sault Ste. Marie STP discharges, extending downstream along the Canadian 
shoreline as far as upper Lake George. 
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TABLE 3.6: MEAN SEDIMENT CHARACTERISTICS ASSOCIATED WITH GUILD B 

DOMINANCE CATEGORIES 



Sediment 
Variable 



Mean* Sediment Characteristic (after transformation) 
Index Range: 0.12 0.12-0.179 0.18-0. 2^*9 0.25 



Water Depth 

Aluminum 

Arsenic 

Cadmium 

Cobalt 

Chromium 

Copper 

Iron 

Mercury 

Manganese 

Nickel 

Lead 

Selenium 

Zinc 

Silt 

Clay 

Silt + Clay 

Coarse (> 1,000 urn) 

Medium (45-1,000 um) 

Fine (< it3 um) 

Dieidrin 

Methoxychlor 

Eh 

Endrin 

Endosulphan II 

Endosulphan 1 

Heptachlor Epoxide 

PC5 

ppDDE 

LOl 

Sol. Extractables 

Total Organic Carbon (TOO 

Total Kjeldahl Nitrogen (TKN) 

PH 

Total Phosphorus 



0.59 


0.73 


0.80 


0.6* 


3.73 


3.92 


3.89 


3.86 


0.*9 


0.7* 


0.81 


1.23 


2.37 


2.69 


2.53 


2.36 


0*72 


0.82 


0.73 


0.51 


1.3« 


l.*9 


1.60 


1.61 


1.15 


1.** 


IM 


1.53 


UM 


t*.32 


*.*1 


*.78 


1.3* 


1.90 


1.83 


2.01 


2.32 


2.*3 


2.52 


3.0* 


QM 


1.12 


1.10 


1.15 


IM 


IM 


t.*6 


LSI 


2.39 


2.9* 


2.77 


2.81 


1.60 


2.03 


2.10 


2.50 


0.62 


0.78 


• 0.66 


0.58 


0.20 


0.28 


0.28 


0.25 


0.6S 


0.90 


0.80 


0.65 


0.09 


0.06 


O.ll 


0.10 


i.O0 


0.77 


0.82 


0.98 


0.5* 


0.78 


0.72 


0.58 


0.88 


0.69 


1.03 


0.7* 


0.61 


0.6a 


0.74 


1.20 


i^78.56 


*57.9* 


*07.25 


262.*3 


0.38 


0.3* 


0.33 


0.67 


0.60 


0.60 


0.65 


0.85 


0.33 


0.33 


0.*7 


0.79 


0.00 


0.05 


0.00 


0.27 


1.35 


1.37 


l.*9 


1.37 


0.Q* 


0.09 


0.18 


o.os 


«>.*7 


6.1* 


6.57 


S.19 


2.70 


3.16 


3.16 


3.26 


1.01 


1.26 


1.3* 


1.91 


2.7* 


2.96 


2.82 


2.7* 


6.27 


5. S3 


6.07 


6.35 


2.52 


2.65 


2.6* 


2.56 



See Table A5.2 for original units and transformation. 
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TABLE 3.7: PREDICTION OF GUILD B DOMINANCE INDEX FROM SEDIMENT 

CHARACTERISTICS 



Sediment 
Variable 



Standard Partial 
Regression Coefficient 



Significance 
Level (p) 



Phosphorus 

Coarse Particulates (> 1,000 urn) 

Eh 

Loss on Ignition (LOl) 

Clay Fraction 

Total Kjeldahl Nitrogen (TKN) 

Heptachlor Epoxide 

Arsenic 

Mercury 

Chromium 

Zinc 

Copper 

Total Organic Carbon (TOO 

Aluminum 

Iron 



-0.177*0 

0.10'*61 
-0.11*636 
-0.8^(502 

0.it7H26 
-0.2547'> 

0.25372 
-0.62972 

0.65302 
-0.30S18 
-0.33140 

0.39308 

1.18693 
-0.i*0969 

0.56982 



0.4106 
0.2675 
0.U21 

Q.OOOl 
0.0026 
0.0954 
0.0079 
0.0196 
0.0002 
C.0396 
0.1594 
0.0909 
C.OOOC 
0.0973 
0.070S 



Source 



df 



Sum of Squares 



Mean Square 



Regression 
Residual 
F = 13.72 
R = 0.88 



14 
55 



0.89083 
0.25501 

R^:s 0.78 



0.06363 

0.00464 

Adjusted R^ = 0.72 
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♦.0 CONCLUSIONS 

The 1985 St. Marys River benthic survey results were not amenable to pooling with those 
of the 19S3 survey. Total densities were found to be generally greater in 1985, moreso at 
some stations than others. Sorters and sorting methods were identified as significant 
contributing factors. Principal components analysis suggested that a faunal shift in 
species associations, either real or sorter-related, had also occurred between the two 
surveys. 

As concluded in 1983, the tubificid community associated with outfall areas was still 
apparent. Nematodes also reached high densities in all but the most severely enriched 
locations. General reductions in pollutant loadings from the Algoma Steel terminal basin 
and the St. Marys Paper operations since 1973 appear to have contributed to minor 
community changes, but were not reflected in major improvements in the poilutional 
status of the benthic community. 

Cluster analysis of sampling stations based on species composition in 1985 showed that 
benthic communities characterized by tolerant species occurred in areas downstrearr, of 
discharges, and other communities characterized by intolerant species occured in areas 
remote from point source effects. Similar impact zones can be derived from a 
community index reflecting dominance of nematodes and immature tubificids without 
setae. Discriminant analysis identified a heavy metal-particle size gradient and a 
pesticide-particle size gradient in sediment quality, which together provide a basis for 
separation of impacted station clusters from each other and from unimpaired stations. 

Impacted stations generally occurred downstream of the Algoma Steel, St. Marys Paper 
and the easterly Sault Ste. Marie sewage treatment plant (STP) discharges. So-r,e 
recovery was apparent with increasing distance downstream from the discharges; 
however, complete recovery was not apparent until the lower section of Lake George. 
The combined impact of effluents from all sources was still evident in Lake George, 
some 2U km downstream from the Algoma and St. Marys Paper discharges. !n general, 
impacted locations were restricted to the Canadian portion of the St. Marys River, *.th 
a clean water fauna characterizing the relatively unindustrialized U.S. shore and all 
portions of the river upstream of pollution sources. 
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The distribution of Hexagenia nymphs was found to be associated with the occurrence of 
visible oil in the substrate of the St. Marys River, Lake Nicolet, and the Lake George 
Channel downstream to the outflow of Lake George. Generally, absence or low densities 
of nymphs coincided with the presence of oil in the sediments. The pattern observed in 
1985 for the distribution of oil and Hexagenia was similar to that reported in 1975, which 
indicated the occurrence of oil as far downstream as the upper half of Lake George. 
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APPENDIX 1 

Historical Assessments of the St. Marys River 

Benthlc Community and Algoma Steel 

Terminal Basin Effluent 
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Figure A1.3 

Distribution and Zones 
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Figure A1.4 

Algoma Steel Terminal Basin 
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Figure A1.5 

Algoma Steel Terminal Basin Effluent 
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Figure A1.6 

Algoma Steel Terminal Basin Effluent 
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Figure A1.7 

Algoma Steel Terminal Basin Effluent 

Sulphides Concentration 
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Figure A1.8 

Algoma Steel Terminal Basin Effluent 
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APPENDIX 2 

Comparability of the 1983 
and 1985 Benthic Surveys 



FIGURE A2.1: SORTER AND TEMPORAL EFFECTS ON LOG TOTAL DENSITY IN 

CORRESPONDING SAMPLES 



3.B-i 



3.*-. 



p,..1 ..F L'"»6 WIIM fTN PV SCOTER 



1. 7 



1 . i75 



Stations 

1 118 

2 SMD 5.0E-2iO 

3 SMD 5,0E-370 
tt SMD 5.0E-520 



S<?rKr 

d = sorter DF 1985 
e = sorter EJ 1985 
m= sorter MS 1983 



66 - 



TABLE A2.1A: DATA FOR EXAMINATION OF SAMPLING METHOD EFFECTS IN 
1985 (no./0.05 m^) 



Total 
Group Station Method Density 



1 1 11,1*56 

1 1 ^,86^* 

J I 3,1*21* 

I i 12,800 

1 2 6,656 

1 2 3,8*0 

1 2 *,28A 

1 . 2 2,560 

1 2 1,792 

J 3 7,509 

J 3 7,765 

J 3 ii,'*iO 

I 3 3.157 

1 3 1,W0 

2 I i,*08 
2 I 1,'*'*0 

2 i 7m 

2 i 2,^96 

2 I l,C2^ 

2 2 752 

2 2 232 

2 .2 800 

2 2 ''96 

2 2 '*00 

2 3 2,059 

2 1 2,731 

2 3 1.760 

2 3 5,909 

2 3 7^1 

I : 1 832 

3 I U656 
3 I 9,66* 
3, 1- 2,002 
3 1 2,*88 
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TABLE A2.1A: 



DATA FOR EXAA^INATION OF SAMPLING METHOD EFFECTS IN 
1985 (no./0.05 m'^) 



Group 



Station 



Method 



Total 
Density 



3 
3 
3 

3 

3 



2 
2 
2 

2 
2 



391 
606 

356 
526 



3 
3 
3 
3 
3 



1,902 
2,177 
2,89^ 
1,3L2 
729 



Station 

1. SMD5.OE-60 

2. SMVl. 0-2130 

3. SMV2.5-1830 



Method 

1. Ponar grab, 200 mm sieve 

2. Ponar grab, 500 mm sieve 

3. Airlift sampler, 200 mm sieve 
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TABLE A2. IB: 



ANALYSIS OF SAMPLING METHOD EFFECTS (DEVICE AND MESH 
SIZE) - LOG DENSITY VALUES (no./0.05 m'^ sample) BY STATION 
AND METHOD (1985 data) 



I. Cell Means 










Station: 


SMD 5.0E-60 


SMU 1.0-2130 


SMU 2.5-1830 


N 




3.65 


IM 


3.05 


15 


Method': 


^0 


P5OO 


A2OO 






3.i»3 


2.93 


3.38 


J5 


SMD 5.0E-60 
SMU 1.0-2130 
SMU 2.5-1830 


3J1 
Ul 
3.36 


3.5^* 
2.69 
2.56 


3.61 
3.33 
3.21 


5 

.5 

5 



Analysis of Variance 



Source of Variation 


Sum of 
Squares 


QL 


Mean 
Square 


F 


Significance 
of F 


Main Effects 

Station 

Method 


5.991 

3.712 
2.279 


2 
2 


L*98 
L856 


20.211 

25.0'»6 

15.376 


0.000* 
0.000* 

Q.OOO* 


Two-Way Interaction 
Station x Method 


0.782 


4 


0.196 


2.639 


0.050** 


Explained 


6.771* 


8 


O.&HJ 


1L'*25 


0.000* 


Residual 


2.668 


36 


0.071* 






TOTAL 


9MI 


1*1* 


0.215 







N = Number of samples. 
* Significant difference at p = 0.05. 
Marginally significant at p = 0.05. 
= Ponar grab, 200 um sieve 



I 



200 
'500 
^200 



= Ponar grab, 500 um sieve 
Airlift sampler, 200 um sieve 
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TABLE A2,1C: 



ANALYSIS OF SAMPLING METHOD EFFECTS (DEVICE ONLY) 
- LOG DENSITY VALUES (no./0.05 m^ sample) BY STATION 
AND METHOD (1985 data) 



1. Cell Means 














Station: 


SMD 5.0E-60 


SMU 


LO- 


2130 


SMU 2.5-1830 


N 




3.71 




3.22 




3.29 


10 


Method^: 


P.10P 
3.i>3 








Ajgo 
3.38 


15 


SMD 5.0E-60 

SMU 1.0-2130 
SMU 2.5-1830 


3.81 

3.II 
3.36 








3.61 
3.3 3 
3.21 


5 
5 



Analysis of Variafw:e 





Sum of 




Mean 




Significance 


Source of Variation 


Squares 


DF 


Square 


F 


of F 


Main Effects 

Station 

Method 


IA25 


3 
2 

1 


0.^*75 
0.706 
O.Qlt* 


5.564 
8.266 

0.161 


0.005* 

0.002* 
0.692 


Two-Way Interaction 
Station x Method 


0.251 


2 


0.125 


1.^69 


0.250 


Explained 


1.676 


5 


3.336 


3.926 


0.010* 


Residual 


2,0^9 


24 


0.085 






TOTAL 


3.725 


29 


0.128 







N = Number of samples. 

* Significant difference at p = 0.05. 
P200 = Ponar grab, 200 urn sieve 
P5OO - P°"3'' S""^*^* ^'^^ '^"^ sieve 



A 



200 



Airlift sampler, 200 um sieve 



ro - 



TABLE A2.2A: DATA FOR^EXAMINATION OF SORTER EFFE 
(no./0.05 m^ sample) IN 1983/85 SAMPLES 


CTSON 


DENSITY 




Total 
Organisms 




Major Taxonomic 


Groups 




Log 


Station 
Sorter 


Chir 


Olig 


Poly 


Nema 


Pele 


Total 

Organisms 


Id 


52 


4 


36 


k 





* 


L72 




56 


6 


m 


2 


m 


4 


L75 




120 


32 


56 


8 


2* 





2.08 


2 4 


k,H9(> 


12 


228 


^^,052 


88 


i 


3.65 




7,888 


m 


512 


6,688 


t&O 


i 


3.90 




5,697 


m 


257 


ii,il6 


256 


112 


3.76 


3d 


7,680 





1,02^ 





6,656 





3.89 




19,968 


1,02't 


9,98'> 





8,70f 





'<.30 




'*,096 





2,0^8 





1,792 





3.61 


4d 


26,62't 


512 


19,968 





6,11*1* 





it.HJ 




1^,336 


512 


12,800 





1,02'f 





i*A6 




13,312 





9,216 





2,560 





f.I2 


le 


91 


J 


*8 


16 





16 


L96 




96 


n 


J& 








S 


L98 




139 


16 


91 


11 


11 


u 


2.14 


2e 


3,390 





320 


2,665 


85 


21 


3.53 




1,015 


9 


tm 


768 


26 


«» 


3.01 




190 


f 


77 





3tt 


26 


2.28 


3e 


256 





236 











2.'»I 




1,82'» 


192 


1,600 





32 





3.26 




IGU 


128 


608 





128 





2.9^* 
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TABLE A2.2A: 



DATA FOR^AMINATION OF SORTER EFFECTS ON DENSITY 
(no./0.05 m sample) IN 1983/85 SAMPLES 





Total 

Organisms 




Major Taxonomic 


GroLps 




Log 


Station 
Sorter 


Chir 


Olig 


Poly 


Nema 


Total 
Pele Organisms 


«e 


3,328 





2,304 





1,024 





3.52 




3,072 





2,560 











3.49 




2,0^*8 





1,536 





512 





3.31 


2m 


'»3.3 


2 


35 








7 


L66 




#0 


2 


21 








11 


r.60 




71.6 


7 


39 








23 


1.85 


3m 


20.3 





19 











1.31 




12 





12 











1.08 




17 





12 











1.23 


Urn 


206.3 


3.3 


203 











2.3! 




132.6 


2.6 


130 











2.12 




238.3 


5.3 


233 











2.38 


Station 
















1 118 




d 


= sorter 


DF 1985 (dissecting microscope 


used) 


2 SMD 5.0E-210 


e 


= sorter EJ 1985 (microscope 


used) 




3 SMD 5.0E-370 


m 


= sorter 


MS 1983 (no microscope used) 




t SMD 3.0E-520 















TABLE A2.2B: 



ANALYSIS OF SORTER EFFECTS ON TOTAL DENSITY IN 1985 
SAMPLES - LOG DENSITY (no./0.05 m^ sample) BY STATION AND 
SORTER 



1. CeU Means 






Station: 




118 


SMD 5.0E-210 






1.94 


3.33 


Sorter ^ 




d 


e 






3.45 


2.82 


lis 

SMD 5.0E- 
5MD 5.0E- 
SMD 5.aE- 


-210 
■370 
-520 


1.85 
3.77 
3.93 
4.24 


2.03 
2.94 
2.87 

3.44 



3.40 



3.84 



N_ 
6 

12 

3 
3 
3 
3 



Analysis of Variance 





Sum of 




Mean 




Significance 


Source of Variation 


Squares 


DF 


Square 


F_ 


of F 


Main Effects 

Station 

Sorter 


14.631 
12.271 
2.360 


4 
3 

i 


3.658 

4.090 
2.360 


36.312 
40.606 
23.431 


0.000* 
O.OOO" 
0.000* 


Two-Way Interaction 
Station x Sorter 


1.368 


3 


0.456 


4.528 


0.018* 


Explained 


15.999 


7 


2.286 


22.690 


o.ooo* 


Residual 


L612 


16 


0.101 






TOTAL 


17.611 


23 


0.766 







N = Number of samiples. 
• Significant difference at p = 0.05. 
^ d ^ sorter DF 1985 
e = sorter EJ 1985 
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TABLE A2.3A; DATA FOR EXAMINATION OF TEMPORAL EFFECTS ON TOTAL 
DENSITY (no./0.03 m sample) 



Total 
Year Station Method Density 



1983 


1 


3 




1 


1 




1 






2 






2 






2 






3 






3 






3 






H 






* 






4 






5 






5 






3 






1 






I 






6 






7 






7 






7 






S 






S 






S 






9 






9 






9 






10 






I# 






m 






,11 






11 






u 




1985 


1 






1 






I 






I 






I 






I 






1 






1 






1 






1 





29 
29 
27 

*3 
39 

71 

20 

12 

17 

204 

133" 

238 

17* 

1*77 

624 

m 

u 

77 

5 

6 

5 

22! 

193 

237 

71 

IQi* 

75 

65 

It 

53 

231 

277 

^♦31 

11,^56 
^,tl6 

12,800 
7,509 
7,765 

3,157 
1,^40 
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TABLE A2.3A: DATA FOR EXAMINATION OF TEMPORAL EFFECTS ON TOTAL 
DENSITY (no./0.05 m^ sample) 



Total 
Year Station Method Density 



1985 


2 


1 








i*,i*96 




2 


I 








7,888 




2 


1 








5,696 




3 


1 








7,680 




3 


I 








25,088 




3 


1 








4,608 




4 


I 








33,280 




'%. 


I 








lf,336 




It 


1 








11,776 




5 


1 








20,224 




5 


1 








18,941^ 




5 


1 








17,408 




6 


1 








3,535 




6 


1 








442 




6 


1 








2,095 




7 


1 








12 




7 


1 








30 




7 


1 








36 




S 


i 








2,502 




& 


1 








8,474 




8 


1 








1,040 




9 


I 








829 




9 


I 








1,640 




9 


1 








9,540 




10 


1 








12,652 




10 


1 








14,410 




10 


1 








13,600 




u 


1 








1,288 




11 


1 








1,408 




11 


1 








640 


Stations 




Me 


thod 








1 SMD 5.0E-60 




1 


Ponar Grab 


, 200 


um sieve 


2 SMD 5.0E-210 




3 


Airii 


ft Sam 


pier, 


200 um sieve 


3 SMD 5.0E-370 














i* SMD 5,0E-520 




Me 


thods 


1 and 


3 are 


shown to be 


5 SMD 5.0E-6'*0 




comparable in 


Table 


A2.1C. 


6 5MU 0.5-200 














7 SMU 0.5-300 














8 SMU 2.5-1300 














9 SMU 2.5-1830 














10 SMU 2.5-2470 














11 SMU 1.0-2130 
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TARLf: A2.1H: ANALYSIS OF TPMPORAl HFrrCJS ON TOTAl nrNSITY. 1981 TO IMJ 

LOC TOTAL DENSITY (no./0.05 oi'' sample) 



I. Cell Mcam 
Station: 

Year: 



SMn 

SMIl 
SMI> 
SMI) 

SMn 

SMIl 
SMI I 
SMIl 
SMH 
SMIl 
SMIl 



5.0E-60 

5.OE-210 

3.0E-J70 

3.0F-320 

).0E-6li0 

0. 5- 200 

0. 5- 300 

2. }- 1 MO 

2.5-1800 

2.3-2^00 

l.0-22f.0 



SMD 5.0F 


SMD 5.0F 


SMH 5. OF 


SMn 5.or 


SMn 5.0F 


SMI) 0.5 


SMU 0.5 


SMIl 2.5 


SMU 2.5 


SMU 2.5 


SMU 1.0 


N 


-60 


-210 


-370 


520 


-MO 


-200 


-100 


-1100 


■ 1800 


-2500 


-2260 




J.19. 


rn 


iy> 


1.26 


IA2 


2.27 


1.05 


2.89 


2.6<t 


2.85 


2.7 5 


6 


1981 


198 5 


N 




















l.7«,. 


3.5<i 


kO 




















IM 


3.71... 






















U*J 


3.77 






















ua 


3.98 






















2.1? 


<l.25 






















ijj 


4.27 






















uu 


1.17 






















0.71. 


l.iT 






















2,Jt 


1.15 


















• 




l.« 


1.37 






















LI6 


o.n 






















2.4» 


J.02 























I 



2. Analysis of Variance 

Sum ol Mean 

Source o( Variation Squar es nF_ l3y??£. 

MainFffetl^ 78.620 II 7.1*7 

SraT7oi^~~ 22.81 » 10 ^-^Xl 

Year 52.70<l I S2.70« 

2-Way Interactions 

Sldlion'iVear'" 8.200 10 0.820 

r«platned 86.820 21 1.1 ll 

Residual 1.S20 51 0.075 

TOTAL 9fl.6<i0 72 1.259 



Significance 

j; oij: 

95.*17 0.000" 

10. 4 58 0.000' 

703.601 0.000* 



10.91(7 
5 5.191 



0.000" 
0.000" 



Number of sjniples. 
N U 
N - 13 
N = 10 
Sif;r>ili' Jill diflerenre at p 



0.05. 



TABLE A2.'»: 1983 BENTHIC SPECIES LOADINGS ON PRINCIPAL COMPONENTS 



Species 



Factor Loadings (Pattern) 
Factor 1 Factor 2 



%of 
Variance 



Hexagenia sp. 
Hyalleia azteca 
Vaivata tricarinata 
Asellus sp. 
Amphipoda sp. indet. 
Bezzia complex 
Mystacides sp. 
Oecetis sp. 
Helisoma sp. 
Protanypus sp. 
Hirudinea 

Rhyacodrilus coccineus 
Parachironomus sp. 
Pagastiella sp. 
Caenis sp. 
Tanytarsus sp. 
Pseudochironomus sp. 
Phryganea sp. 
Dicrotendipes sp. 
Pisidium sp. 
Amnicola limosa 
Paramerina sp. 
Polycentropus sp. 
Lumbricidae sp. indet. 
Epicocladius sp. 
Spirosperma ferox 
Procladius sp. 
Demicryptochirononnus sp. 
Phyiocentropus sp. 
Physa gyrina 
Sphaerium sp. 
Ischnura sp. 
Phaenopsectra sp. 
Microtendipes sp. 
Polypedilunn sp. 
Sialis sp. 
Molanna sp. 
Enallagma sp. 
Pianaria sp. 
Pisidium casteranum 
Sphaerium simile 
Pisidium compressum 
Thienemannimyia complex 
Tubifex harmoni 



0.984 
0.956 
0.951 
0.936 
0.935 
0.918 
0.918 
0.918 
0.918 
0.918 
0.918 
0.918 
0.918 
0.918 
0.884 
0.884 
0.884 
0.884 
0.868 
0.851 
0.798 
0.785 
0.773 
0.692 
0.692 
0.627 
0.623 
0.612 
0.269 
0.269 
0,269 
0.269 
0.436 
0.643 
0.678 
0.000 
0.000 
0.000 
0.000 
-0.308 
0.000 
0.000 
-0.326 
0.000 



0.000 
0.000 
0.000 
0.349 
0.000 
-0.340 
-0.340 
-0.340 
-0.340 
-0.340 
-0.340 
-0.340 
-0.340 
-0.340 
0.453 
0.453 
0.453 
0.453 
0.000 
-0.354 
0.274 
-0.371 
0.300 
-0.360 
-0.360 
-0.328 
0.000 
-0.355 
0.948 
0.948 
0.948 
0.948 
0.888 
0.757 
0.726 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 



0.968 
0.914 
0.904 
0.998 

0.874 

0.958 

0.958 

0.958 

0.958 

0.958- - 

0.958 

0.958 

0.958 

0.958 

0.987 

0.987 

0,987 

0.987 

0.753 

0.850 

0.712 

0.754 

0.688 

0.608 

0.608 

0.501 

0.388 

0.501 

0.971 

0.971 

0.971 

0.971 

0.979 

0.986 

0.987 

0.000 

0,000 

0.000 

0.000 

0.095 

0.000 

0.000 

0.106 

0,000 
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TABLE A2.*: 



1983 BENTHIC SPECIES LOADINGS ON PRINCIPAL COMPONENTS 





Factor Loadings 


(Pattern) 


%of 
Variance 


Species 


Factor 1 


Factor 2 


Limnodrtlus udekemianus 


0.000 


0.000 


0.000 


Immature Tubificidae w/o hair setae 


0.000 


0.000 


0.000 


Potamothrix vejdovskyi 


0.000 


0.000 


0.000 


Simile striatium 


0.000 


0.000 


0.000 


Enchytraeidae 


0.000 


0.000 


0.000 


Proioma rub rum 


0.000 


0.000 


0.000 


Limnodrilus hoffmeisteri 


-0.279 


0.000 


0.078 


Aulodriius americanus 


0.000 


-0.321 


0.103 


Larsia sp. 


0.329 


-0.1^56 


0.316 


Stylodrilus heringianus 


0.000 


0.000 


0.000" 


Cricotopus sp. 


0.000 


0.000 


0.000 


Tanypus sp. 


0.000 


0.000 


0.000 


Aulodriius piqueti 


0.000 


0.000 


0.000 


Immature Tubificidae with hair setae 


-0.329 


0.000 


0.108 


Vejdovskyella comata 


0.000 


0.000 


0.000 


Ablabesmyia sp. 


0.000 


0.000 


0.000 


Quistadrilus multJsetosus 


0.000 


0.000 


0.000 


Glossoscolecidae 


0.000 


0.000 


0.000 


Uncinais uncinata 


0.000 


0.000 


0.000 


Stylaria lacustris 


0.000 


0.000 


0.000 


Specaria josinae 


0.000 


0.000 


o.oco 


Chironomus sp. 


0.000 


0.000 


0.000 


Cladopelma sp. 


0.000 


0.000 


0.000 


Cryptochironomus sp. 


0.393 


.O.ti60 


0.366 


Cryptotendipes sp. 


0.000 


0.000 


0.000 


Manayunkia speciosa 


0.000 


0.000 


0.000 


Slavina appendiculata 


0.000 


0.000 


0.000 


Glyptotendipes sp. 


0.000 


0.000 


0.000 


Harnischia sp. 


0.000 


0.000 


0.000 


Nais pardalis 


0.000 


0.000 


0.000 


Lumbriculus variegatus 


0.000 


0.000 


0.000 


Gammarus sp. 


0,000 


0.000 


0.000 


Paraiauterborniella sp. 


0.000 


0.000 


0.000 


Nais bretscheri 


o.oco 


0.000 


o.ooc 


N. variabilis 


0.000 


0.000 


0.000 


Lymnaea sp. 


0.000 


0.000 


0.000 


Rheotanytarsus sp. 


0.000 


0.000 


0.000 


Stempellina sp. 


0.000 


0.000 


0.000 


Stictochironomus sp. 


0.000 


0.000 


0.000 


Ripistes parasitica 


0.000 


0.000 


0.000 


Isochaetides curvisetosus 


0.000 


0.000 


0.000 


Piscicola punctata 


0.000 


0.000 


0.000 


Heterotrissocladius sp. 


0.000 


0.000 


0.000 


Orthocladius sp. 


0.000 


0.000 


0.000 
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TABLE A2M 1983 BENTHIC SPECIES LOADINGS ON PRINCIPAL COMPONENTS 



Species 



Psectrocladius sp. 
Thienemanniella sp. 
Diamesa sp. 
Monodiamesa sp. 
Crangonyx sp. 
Collembola sp. indet. 
Potthastia sp. 
Ephemera sp. 
Hydra sp. 

Paraleptophlebia sp. 
Valvata sincera 
Stenonema sp. 
Ephemeroptera sp. indet. 
Pristina osborni 
Nais barbata 
Piquetiella sp. 
Agraylea sp. 
Oxyethira sp. 
Ceraclea sp. 
Lirceus sp. 
Nectopsyche sp. 
Nematoda sp. indet. 
Pristina foreli 
Ptilostomis sp. 
Stylaria flossularis 
Campeloma decisum 
Aulodrilus pluriseta 
Hydracarina sp. indet. 
Baetis sp. 



Factor Loadings 


(Pattern) 


%oi 


Factor 1 


Factor 2 


Variance 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


o.ooo 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


o.oco 


0,000 


0.000 


0.000 


0.000 


0.000 


o.oov 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0,000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 



% of Variance 



23.36if 



9.395 



See plot of species loadings on Factors 1 and 2 in Figure A2.2. 
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TABLE A2.3: 1983 BENTHIC SPECI 


ES LOADINGS ON 


PRINCIPAL C 


:OMPONENT5 




Factor Loadings 


(Pattern) 


%of 


Species 


Factor 1 


Factor 2 


Variance 


Vaivata sincera 


0.985 


0.000 


0.970 


Pisidium sp. 


0.975 


0.000 


0.951 


Tany tarsus sp. 


0.975 


0.000 


0.951 


Pagastiella sp. 


0.969 


0.000 


0.939 


Bezzia complex 


0.96^ 


0.000 


0.929 


Microtendipes sp. 


0.961 


0.000 


0.924 


Polycentropus sp. 


0.957 


0.000 


0.916 


Cryptotendipes sp. 


0.953 


0.000 


0.908 


Hexagenia sp. 


0.952 


0.000 


0,906 


Paraiauterborniella sp. 


0.951 


0.000 


0.90V 


Glyptotendipes sp. 


0.9if8 


0.000 


0.899 


Caen is sp. 


0.948 


0.000 


0.899 


Nectopsyche sp. 


0.9*S 


0.000 


0.899 


Sialis sp. 


0.9^*8 


0.000 


0.899 


Campeloma decisum 


0.944 


0.000 


0.891 


Mystacides sp. 


0.936 


0.000 


0.876 


Aulodrilus pluriseta 


0.929 


0.000 


0.863 


Psectrocladius sp. 


0.843 


0.000 


0.711 


Harnischia sp. 


0.834 


0.344 


0.814 


Amnicoia limosa 


0.796 


0.000 


0.634 


Stempeilina sp. 


0.785 


0.459 


0.827 


Aulodrilus americanus 


0.678 


0.510 


0.720 


Crictopus sp. 


0,674 


0.000 


0.454 


Lirceus sp. 


0.593 


0.447 


0.551 


Lymnaea sp. 


0.000 


0.948 


0.899 


Monodiamesa sp. 


0.000 


0.934 


0.872 


Ephemera sp. 


0.000 


0.889 


0.790 


Rheotanytarsus sp. 


0.353 


0.701 


0.616 


Oecetis sp. 


0.000 


0.700 


0.490 


Specaria sp. 


0.613 


0.629 


0.771 


Nematoda sp. indet. 


0.000 


-0.614 


0.377 


Immature Tubificidae w/o hair setae 


0.000 


-0.574 


0.329 


Baetis sp. 


0.000 


0.625 


0.391 


Campeloma decisum 


0.000 


0.625 


0.391 


Styiaria fossularis 


0.000 


0.625 


0.391 


Ceraclea sp. 


0.000 


0.625 


0.391 


Nais bretscheri 


0.000 


0.625 


0.391 


Thienemannieila sp. 


0.000 


0.625 


0.391 


Pristina foreli 


0.000 


0.625 


0.39 J 


Stenonema sp. 


0.000 


0.625 


0.391 


Isochaetides cirvisetosus 


0.000 


0.625 


0.391 


Pique ti el la sp. 


0.000 


0.625 


0.391 



- 80 - 



TABLE A2.5; 1983 BENTHIC SPECIES LOADINGS ON PRINCIPAL COMPONENTS 



Species 



Factor Loadings (Pattern) 
Factor 1 Factor 2 



Gammarus sp. 
Ephemeroptera sp. indet. 
Stictochironomus sp. 
Hyallela azteca 
Parachironomus sp. 
Uncinais uncinata 
Paraleptophlebia sp. 
Pseudochironomus sp. 
Stylaria lacustris 
Epoicocladius sp. 
Hydra sp. 

Heterotrissocladius sp. 
Crangonyx sp. 
Ptilostomis sp. 
Helisoma sp. 
Hydracarina sp. indet. 
Limnodriius udekemianus 
Quistadrilus multisetosus 
Agraylea sp. 
Abiabesmyia sp. 
Linnnodriius hoffmeisteri 
Aseilus sp. 
Polpedilum sp. 
DicrotendiF>es sp. 
Nais variabilis 
Larsia sp. 
Piscicoia punctata 
Potthastia sp. 
Physa gyrina 
Oxyethira sp. 
Enchytraeidae sp. indet. 
Sphaerium sp. 
Phaenopsectra sp. 
G losses CO lee idae sp. indet. 
SJavina appendiculata 
Orthocladius sp. 
Cladopelma sp. 
Cryptochironomus sp. 
Procladius sp. 
Auiodrilus piqueti 
Thienemannimyia sp. 
Demicryptochironomus sp. 



0.000 
0.000 
0.000 
0.000 

0.000 
0.000 
0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

O.SQM 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.6S5 

0.000 

Q.5U0 

0.555 

0.000 

0.000 

0.000 

0.330 

0.*^8 

0.000 

0.3S9 

0.000 



0.625 
0.625 
0.625 
0.632 
0.56f 
0.575 
0.575 
0.575 
0.575 
0.569 
0.595 
0.629 
0.515 
-0.388 
-0.'f28 
-0.'f28 
-O.'flS 
-0.^*18 
-0.^18 
-0.^76 
-0.^*50 
-0.282 
0.000 
-0.^13 
-0.50t 
0.000 
0.000 
0.000 
0.'»06 
0.ttl3 
0.000 
0.620 
0.507 
0.000 
-0.279 
0.000 
0.000 
0.000 
0.000 
0.329 
0.286 
0.000 



%of 


Variance 


0.391 


0.391 


0.391 


0.399 


0.318 


0.331 


0.331 


0.331 


0.331 


0.32T 


0.351^ 


0.396 


0.265 


0.151 


0.183 


0.183 


0.175 


0.175 


0.175 


0.227 


0.203 


0.080 


0.25^* 


0.171 


0.25'f 


0.000 


0.000 


0.000 


0.165 


0.171 


0.^*69 


0.382* 


0.5^*9 


0.308 


0,078 


0.000 


0.000 


0.109 


0.201 


0.108 


0.233 


0.000 
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TABLE A2.5: 1983 BENTHIC SPECIES LOADINGS ON PRINCIPAL COMPONENTS 





Factor Loadings 


(Pattern) 










%of 


Species 


Factor I 


Factor 2 


Variance 


Manayunkia speciosa 


0.000 


0.000 


0.000 


Ripistes parasitica 


0.000 


0.000 


0.000 


Pristina osborni 


0.000 


0.000 


0.000 


Stylodrilus heringianus 


0.000 


0.000 


0.000 


Vejdovskyella comata 


0.000 


0.000 


0.000 


Nais barbata 


0.000 


0.000 


0.000 


N. pardalis 


0.000 


0,000 


0.000 


Diamesa sp. 


0.^66 


-0.31^ 


0.316 


Spirosperma ferox 


0.000 


-0.4^5 


0.198 


Prostoma rubrum 


0.265 


-0.36 if 


0.203- 


Chironomus sp. 


0.000 


0.000 


0.000 


Immature Tubificidae with hair setae 


0.000 


-0.388 


0.151 


Collembola sp. indet. 


0.000 


0.000 


0.000 


Lumbricilus variegatus 


0.000 


0,000 


0,000 


Pisidium casertanum 


0.000 


0,000 


0.000 


Protanypus sp. 


0.000 


0.000 


0.000 


Sphaerium simile 


0.000 


0.000 


0.000 


Enallagma sp. 


0.000 


0.000 


0.000 


Ishnura sp. 


0.000 


0.000 


0.000 


^haerium striatium 


0.000 


0.000 


0.000 


Tanypus sp. 


0.000 


0.000 


0.000 


Pisidium compressum 


0,000 


0.000 


0.000 


Paramerina sp. 


0.000 


0.000 


0.000 


Valvata tricarinata 


0.000 


0.000 


0.000 


Rhyacodrilus coccineus 


0.000 


0.000 


0.000 


Lumbricidae sp. indet. 


0.000 


0.000 


0.000 


Molanna sp. 


0.000 


0.000 


0.000 


Tubifex harmani 


0.000 


0.000 


0.000 


Phryganea sp. 


0.000 


0.000 


0.000 


Phylocentropus sp. 


0.000 


0,000 


0.000 


Hirudinea sp. indet. 


0.000 


0.000 


0.000 


Amphipoda sp. indet. 


0.000 


0.000 


0.000 


Potamothrix vejdovskyi 


0.000 


0.000 


0.000 



% of Variance 



23.196 



18.632 



See plot of species loadings on Factors 1 and 2 on Figure A2.2. 
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FIGURE A2.2: PLOTS OF SPECIES LOADINGS ON PRINCIPAL COMPONENT FACTORS: 
A - 1983, B - 1985, * - Single Species, $ ^ Multiple Species 



APPENDIX 3 
Ouster Analysis of I9i5 Benthic Data 
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TABLE A3.1: CONCORDANCE OF CLUSTER ANALYSIS SOLUTIONS USING 

DIFFERENT CLUSTERING METHODS 



WARD: Hierarchical .Agglomeration - Ward's Method 

WAVG: Hierarchical Agglomeration - Within Group Average Link 

BAVC: Hierarchical Agglomeration - Between Group Average Link 

KMEANS(7): K-Means Partition - 7 Cluster Solution 

KMEANS(S): K-Means Partition - S Cluster Solution 

GHOW: Hierarchical Polythetic Division Howard-Harris Method 
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